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Abstract 

In three-dimensional O(A^) models, we investigate the low-momentum be- 
havior of the two-point Green's function G{x) in the critical region of the 
symmetric phase. We consider physical systems whose criticality is charac- 
terized by a rotational-invariant fixed point. Several approaches are exploited, 
such as strong-coupling expansion of lattice non- linear O(A^) a models, 1/N- 
expansion, field-theoretical methods within the (j)'^ continuum formulation. 

Non-Gaussian corrections to the universal low-momentum behavior of 
G{x) are evaluated, and found to be very small. 

In non-rotational invariant physical systems with 0(A^)-invariant interac- 
tions, the vanishing of space-anisotropy approaching the rotational-invariant 
fixed point is described by a critical exponent /3, which is universal and is 
related to the leading irrelevant operator breaking rotational invariance. At 
N = oo one finds p = 2. We show that, for all values of > 0, p ~ 2. 

PACS numbers: 05.70.Jk, 64.60.Fr, 75.10.Hk, 75.40.Cx 
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I. INTRODUCTION 



Three-dimensional O(A^) models describe many important critical phenomena in nature. 
The statistical properties of ferromagnetic materials are described by the case N = 3, where 
the lagrangian field represents the magnetization. The helium superfluid transition, whose 
order parameter is the complex quantum amplitude of helium atoms, corresponds to = 2. 
The case = 1 (i.e. Ising-like systems) describes liquid-vapour transitions in classical fluids 
or critical binary fluids, where the order parameter is the density. O(A^) models in the limit 
^ describe the statistical properties of long polymers. 

The critical behavior of the two-point correlation function G{x) of the order parameter 
is relevant in the description of critical scattering observed in many esperiments, typically 
neutron scattering in ferromagnetic materials, light and X rays in liquid-gas systems, ... In 
Born's approximation the cross section Tfi for incoming particles (i.e. neutrons or photons) 
of momentum pi and final outgoing momentum P/ is proportional to the component k = 
Pf — Pi of the Fourier transform of G{x) 

Tfi(xG{k = pf -pi). (1) 

As a consequence of the critical behavior of the two-point function G{x) at T,,, 

- (2) 

the cross section for A; ^ (forward scattering) diverges as T ^ T^. When strictly at 
criticality the relation (0) holds at all momentum scales. In the vicinity of the critical 
point where the relevant correlation length ^ is large but finite, the behavior @) occurs 
for A ^ A; ^ 1/^, where A is a generic cut-off related to the microscopic structure of the 
statistical system, for example the inverse lattice spacing in the case of lattice models. At 
low momentum, -C 1/^, experiments show that G{x) is well approximated by a Gaussian 
behavior, 

where Mq ~ 1/^ is a mass scale defined at zero momentum (for a general discussion see e.g. 
Ref. i). 

We will specifically consider systems with an 0(A^)-invariant Hamiltonian in the sym- 
metric phase, where the 0{N) symmetry is unbroken. Furthermore, we will only consider 
systems with a rotationally-symmetric fixed point. Interesting members of this class are 
systems defined on highly symmetric lattices, i.e. Bravais or two-point basis lattices with a 
tetrahedral or larger discrete rotational symmetry. 

In this paper we focus on the low-momentum behavior of the two-point correlation func- 
tion of the order parameter, which coincides with the lagrangian field, in three-dimensional 
0{N) models. We want to estimate the deviations from Eq. (^) in the critical region of 
the symmetric phase, i.e. for < T/Tc — 1 <^ 1, and in the low-momentum regime, i.e. 
fc^ ^ Mq. We focus on two quite different sources of deviations: 
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(i) Scaling corrections to Eq. (^), depending on the ratio k'^/MQ, and reflecting the 
non-Gaussian nature of the fixed point. 

(ii) Non-rotationally invariant scaling violations, reflecting a microscopic anisotropy in 
the space distribution of the spins. This phenomenon may be relevant, for example, in the 
study of ferromagnetic materials, where the atoms lie on the sites of a lattice giving rise to a 
space anisotropy which may be observed in neutron-scattering experiments. In these systems 
the anisotropy vanishes in the critical limit, and G{x) approaches a rotationally invariant 
form. It should be noticed that this phenomenon is different from the breakdown of the 
O(A^) symmetry in the interaction, which has been widely considered in the literature [Q. 

In our study of the critical behavior of the two-point function of the order parameter 
G{x) we will consider several approaches. We analyze the strong-coupling expansion of G{x), 

G{x)^{s{x)-m), (4) 
for the lattice O(A^) non-linear a model with nearest-neighbor interaction 

Sl = -N(3 Y1 ^--^y^ (5) 

links {xy) 

which we have calculated up to 15th order on the simple cubic lattice and 21st order on the 
diamond lattice. We also perform a detailed study using the 1/A^-expansion, whose results, 
beside clarifying physical mechanisms, are also useful as benchmarks for the strong-coupling 
analysis. Moreover we compute the first few non-trivial terms of the e-expansion and of 
the (7-expansion (i.e. expansion in the four-point renormalized coupling at fixed dimension 
c? = 3) of the two-point function for the corresponding (p^ continuum formulation of O(A^) 
models: 

We recall that non-linear a models and 0^ models possessing the same internal symmetry 
O(A^) describe the same critical behavior. By universality our study provides information 
on the behavior of the physical systems mentioned above in the critical region of the high- 
temperature phase. A short report of our study can be found in Ref. 0. 

The first systematic study of the critical behavior of G{x) is due to Fisher and 
Aharony They computed G{x) in the e-expansion up to terms O(e^) and in the 

large- expansion to order 1/A^ [Q; moreover some estimates of the non-Gaussian correc- 
tions for A^ = 1 and A^ = 3 were derived from strong-coupling series for various lattices . 
Their calculations confirmed experimental observations that non-Gaussian corrections are 
small in the low-momentum region. 

In this paper we reconsider the problem of determining the non-Gaussian correction to 
G{x) in the low-momentum regime using the different approaches we mentioned above. Our 
strong-coupling analysis of the A^-vector model on the simple cubic and diamond lattice 
leads to rather accurate results which considerably improve earlier calculations. This is 
achieved essentially for two reasons: longer strong-coupling series are available, and, more 
importantly, we consider improved estimators which allow more stable extrapolations to the 
critical limit. Our strong-coupling analysis is integrated and supported by results obtained 
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from the 0{1/N), e- and (^-expansion. We compute the expansion of G{x) up to four-loops, 
0{g'^), in fixed dimension c? = 3 and we extend the results of 0] by calculating the next three- 
loop term O(e^). The results of the various approaches are reasonably consistent among each 
other: the ^f-expansion and the analysis of the strong-coupling series provide in general the 
most precise estimates, together with the expansion for ^ 16. The e-expansion is 
somewhat worse but still consistent, perhaps because of the limited number of terms (one 
term less than in the (/-expansion). 

We also discuss the space anisotropy in G{x) induced by the lattice structure. For the 
class of systems we consider, G{x) becomes rotationally-invariant at criticality: when j3 — *• 
(3ci so that Mg — * 0, the anisotropic deviations vanish as M^, where p is a universal critical 
exponent. From a field-theoretical point of view, space anisotropy is due to non-rotationally 
invariant 0(A^)-symmetric irrelevant operators in the effective Hamiltonian, whose presence 
depends essentially on the symmetries of the physical system or of the lattice formulation. 
The exponent p is related to the critical effective dimension of the leading irrelevant operator 
breaking rotational invariance. On (i-dimensional lattices with cubic symmetry the leading 
operator has canonical dimension d+2. In the large- limit, where the canonical dimensions 
determine the scaling properties, one finds p = 2 with very small 0{1/N) corrections. A 
strong-coupling analysis supported by a two-loop e-expansion and three-loop (/-expansion 
computation indicate that p remains close to its canonical value for all iV > 0, with deviations 
of approximately 1% for small values of A^. It should be noted that the exponent p which 
controls the recovery of rotational invariance is different from cu, the leading subleading 
exponent, since they are related to different irrelevant operators. This means — and this 
may be of relevance for numerical calculations — that the recovery of rotational invariance 
is unrelated to the disappearance of the subleading corrections controlled by tu: in practice, 
as p ~ 2 while 0.5 ^ ^ 1 [^0, rotational invariance is recovered long before the scaling 
region. 

We also investigated the recovery of rotational invariance in two-dimensional models. On 
the square lattice, for iV = 1 (Ising model) and > 3, we show that p = 2. This leads us 
to conjecture that p = 2 holds exactly for all two-dimensional models on the square lattice. 
Similarly we conjecture that p = 4 (resp. p = 3) are the exact values of the exponents for 
the triangular (resp. honeycomb) lattice. A Monte Carlo and exact-enumeration study |§] 
for A^ = on the square lattice is consistent with this conjecture. 

We should mention that our results on space anisotropy are also relevant in the discussion 
of the linear response of the system in presence of an external (anisotropic) field. 

The paper is organized as follows: 

In Sec. ^ we fix the notation and introduce a general parametrization of G{x) which 
includes the off-critical and non-spherical dependence. 

In Sec. |Tl| we analyze the critical behavior of G{x) at low momentum. We present 
calculations based on various approaches: 1/A^-expansion (up to 0(1/A^)), (/-expansion (up 
to 0{g'^)), e-expansion (up to 0(e^)), and an analysis of the strong-coupling expansion of 
G{x) on the cubic and diamond lattice. 

In Sec. 13 the anisotropy of G{x) is studied in the critical region. We present large- A?" 
and 0(1/A^) calculations on various lattices, and a strong-coupling analysis of some non- 
spherical moments of G{x) on cubic and diamond lattice. Again, the analysis of the first 
non-trivial terms of the ^f-expansion and the e-expansion is presented. Anisotropy in G{x) 
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is also studied in two-dimensional O(A^) models. 

In App. ^ we collect explicit formulae for the large-A^ limit of G{x) for the nearest- 
neighbor A^-vector model on the cubic, f.c.c. and diamond lattices. 

In Apps. and we present some details of our 0{1/N) calculations. 

In App. 0we present the 15th-order strong-coupling expansion of the two-point function 
on the cubic lattice. 

In App. ^ we report the 21st-order strong-coupling series of the magnetic susceptibility 
and of the second moment of G{x) on the diamond lattice for A^ = 1, 2, 3. 



II. THE TWO-POINT GREEN'S FUNCTION 

A. Hypercubic lattices 

In this Section we discuss the general behaviour of the two-point spin-spin correlation 
function in lattice O(A^) non- linear a models. We consider a generic Hamiltonian defined 
on a hypercubic lattice 

n = - ^J{x - y)s^ ■ Sy (7) 
x,y 

where the sum runs over all lattice sites. We will later extend our analysis to other lattices. 
Let us define 

li^ik) = 2 [j{k) - J(0)] (8) 

where J{k) denotes the Fourier transform of J{x). In spite of the notation, we are not 
assuming that k (k) is a sum of the type J2iif{k^). We consider models for which, by a 
suitable normalization of the inverse temperature (3, one finds 

k\k) = e + o{k^), (9) 

so that the critical limit is rotationally invariant. Moreover we make the following assump- 
tions: 

2 

(1) The interaction J(x) is short-ranged so that k is continuous; 

2 

(2) The function J(x) (and thus also k ) is invariant under all the symmetries of the 
lattice; 

2 

(3) The interaction is ferromagnetic, so that A; = only for A; = in the Brillouin zone. 
Beside the leading (universal) rotationally-invariant critical behaviour, we are interested 

in understanding the effects of the lattice structure on the two-point function and the re- 
covery of rotational invariance. For this reason, our analysis must take into account the 

— 2 

irrelevant operators which break rotational invariance. It is natural to expand k {k) in 
multipoles by writing 

oo Pi 

k\k) = Y.Y.e^i\k')Q^^\k). (10) 

1=0 p=l 
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Here the functions Q2i\k) are multipole combinations which are invariant under the sym- 
metries of the lattice. Their expressions can be obtained from the fully symmetric traceless 
tensors of rank 21, T^^" "^'(/c) 0,0, by considering all the cubic-invariant combinations, 
which can be obtained by setting equal an even number of indices larger than or equal to 
four and then summing over them. Odd-rank terms are absent in the expansion (0) because 
of the parity symmetry x — > —x. Moreover, there is no rank-two term, i.e. Q2{k) = 0, due to 
the discrete rotational symmetry of the lattice. The summation over p in Eq. (p!0| ) is due to 
the fact that, for given /, there are in general many multipole combinations which are cubic 
invariant. The number pi depends on the dimensionality d of the space; more precisely it is 
a non-decreasing function of d. It can be computed from the following generating function 
(for the derivation see Appendix |D|) 



d 



Hvit' = n Y-^- (11) 

Z=0 m=2 ^ ^ 

For notational simplicity, we will suppress the explicit dependence on p in all the following 
formulae, but the reader should remember that it is understood in the notation. 

Let us give the explicit expressions of Q2i\k) for the first few values of /. We set 
Qo{k) = 1. For / = 2 there is only one invariant combination, i.e. p2 = 1, which can be 
derived from 

Tf"'\k) = k'^k^k-'k^ - — ^Perm„.^5 f^^^FA;^ 

[d + A) ^ 

where Permc^.. .„„(...) represents the sum of the non-trivial permutations of its arguments. 
One then defines 

Q^{k) = y: rr' = k'- j^ikr, (13) 

where the notation fc" = k^ is used. For / = 3, = 1 for all d > 2. From 

yM-a/37<5^^) = k^k^k^k^k^k^ - — ^Perm,/3^5^, (6'"'k''kf'k''k^) + (14) 

(tt + 8) ^ ^ 



one finds 



In (i = 2 it is easy to verify that Qelk) = so that p^ = 0. For I = 4 and d > 3 two different 
Qg'\k) can be extracted from the corresponding tensor j'^i --"*; 
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Q« (k) = E Tg™™, (16) 
Qf\k) = J2T^'""""""'. (17) 

When d = 2,3 the two combinations are not independent. Indeed Qg^^ = 2Q^^^ so that 
P4 = 1. Higher values of / can be dealt with similarly. 

In order to study the formal continuum limit of the Hamiltonian defined in (|^), we 
expand e2i{k'^) in powers of k"^. We write (the sum over different multipoles with the same 
value of / being understood in the notation) 

oo oo 

k\k) = J2Y. e2i,mik^rQ2iik), (18) 

1=0 m=0 



where eo,o = and eo,i = 1. Inserting back in Eq. (j^) one sees that Eq. (|T8D represents an 
expansion in terms of the irrelevant operators 

02iA^) = s{x) ■ D"'Q2i{d)s{x), (19) 

where □ = J^^^d^- The leading operator that breaks rotational invariance is the four- 
derivative term 

O^ix) = 04,o(a;) = s{x) ■ Qi{d)s{x), (20) 

which has canonical dimensions d + 2. 

Let us now consider the Green's function 

G{x;P) = {so-s,), (21) 

and its Fourier transform G{k; j3). We define a zero-momentum mass-scale Mg{(3) by 

MaW) . ^ (22) 

where Cg(/5) is the second-moment correlation length 

^^^^^ - Yd E.G(x;/5) • (''^ 

Since there is a one-to-one correspondence between Mg{i3) and (3, one may consider G{k] j3) 
as a function of Mg{P) instead of /?. Indeed, for the purpose of studying the critical limit, 
it is natural to consider G{k; P) as a function of k and Mq. In complete analogy to our 
discussion of J{k), we analyze the behaviour of G{k,MG) in terms of multipoles (again 
a sum over different multipole combinations with the same value of / is understood, see 
Eq. (0)): 

oo 

G~\k,MG) = Y.g2iiy,MG)Q2i{k), (24) 

1=0 
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where y = k"^ /Mq. Notice that 

Q2i{k) = Q2i{k/MG)M^l (25) 

For the purpose of studying the universal properties of the critical limit of in 
which Mg — >■ keeping k/Mc fixed, it is important to understand the behavior of the 
functions g2i{y,MG) when Mq 0. The naive limit does not exist. However, as long 
as the contributions to Mq) are originated by the insertion of individual (irrelevant) 

operators without any mixing among different operators with the same symmetry properties, 
one can apply standard results in renormalization theory. In this case one can establish some 
universal properties. For a generic choice of J{k) this holds only for the functions go{y, Mq) 
and g^lyjMc)- Indeed for higher values of / there are mixings among different operators 
which make the renormalization of the functions g2i{y, Mq) more complicated. Consider, for 
instance, the case / = 3 in the large- limit, where the operators have canonical dimensions. 
In this case terms proportional to Q^i^k) are depressed as Mq, while terms proportional to 
Q^ik) are depressed as Mq. However it is easy to see that the multipole decomposition of 
Qi{kY, which is also depressed as Mq, contains a term of the form k'^QQ{k). This means 
that there are two operators contributing to gQ^kjMc), 04^0(2;)^ and Oq^^x). An analogous 
argument applies to higher values of /. Notice that for the particular case of / = 3 the mixing 
should disappear in the limit y — > 0: thus for Mq — > ge{0, Mc) can be directly related to 
the renormalization properties of the operator OQfl{x). 

For / = and / = 2 standard results of renormalization theory show that, if Z2i{Mg) = 
(72/(0, Mc), the following limit exists: 

g2i{y,MG) ^ . . , . 

where g2i{y) is a smooth function which is normalized so that ^2/(0) — 1- The function g2i{y) 
is universal, in the sense that it is independent of the specific Hamiltonian. 

The function g^ly) can also be obtained by considering the linear response of the system 
to an external field possessing the corresponding symmetry properties. One considers the 
one-particle irreducible two-point function with an insertion of a 02i,o{x) operator at zero 
momentum, i.e. 

ro,(a:, Mg) = J dz {02i,o{z) s(0) ■ six))'"-' (27) 
and the corresponding Fourier transform ro2i{k, Mq)- Setting 



the following limit exists 



To,Xk,MG) ^ ^ ( )g (^)_ (29) 



For / = 2 the function defined by the previous equation coincides with that defined in 
(B^); moreover for Mg — > 0, Z4^{Mg)/ Zi{MG) is a finite (non-universal) constant, meaning 
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that both quantities have the same singular behavior for Mq — ^ 0. For higher values of 
/, formula ( p9|) still holds, but there is no easy relation between g2i{y) and g2i{y,MG) as 
defined in Eq. (|2^), at least for generic Hamiltonians. Indeed, at least in principle, one may 
consider specific forms of J{k) enjoying the property that all contributions g2n{k-, Mg) with 
< n <l vanish in the critical limit, for a given value of I. In lattice quantum field theory 



this is essentially the spirit of Symanzik's improvement program |Tl||. In this case formula 
(PB|) is valid for / = I and the corresponding function g2i{y) coincides with that defined by 
Eq. m. 



The functions g2i{y) defined in Eq. ( P^ ) have a regular expansion in y around y = 0: 

g2l{y) = 1 + C2Z,l2/ + C2/,2l/^ + . . . . (30) 

Co,i = 1 due to the definition of the second-moment correlation length. 

The renormalization constant Z2i{Mg) is instead non-universal. For Mq — > it behaves 

as 

Z2i{MG)^Z2iMa''''\ (31) 

where ri2i is a critical exponent which depends only on the spin of the representation (i.e. 
it does not depend on the additional index p which has always been understood in the 
notation, see Eq. ([T0|)), and 121 is a non-universal constant which depends on the lattice 
and on the Hamiltonian (and the additional index p). An analogous expression is valid for 
Zi^Mc) (and for Z^i for the special Hamiltonians we have discussed before): for Mq 
we have Z4^{Mg) ~ Z4^Mq^^. For / = 0, as a consequence of our definitions, 

Zo{Mg) ~ M^-^ (32) 

where rj is the standard anomalous dimension of the field. More generally a2i = rj — rj2i is 
the anomalous dimension of the irrelevant operator 02ifi{x). 

In two dimensions the renormalization constants diverge only logarithmically and thus 
we write for I 7^ 



Z2i{Mg)^Z2i{\^Mg: 



72! 



1 + 



1 



InM, 



G 



(33) 



The anomalous dimensions 72^ are universal while the prefactor 121 depends on the details 
of the interaction. 

We can now discuss the critical limit of Eq. (|2^). Using the previous formulae we can 
write for Mq 

^ Uy) + "rot. inv. sublead." + ^M^+''-''-^4(?/)Q4(A:/Mg) + . . . (34) 

where "rot. inv. sublead." indicates rotationally-invariant sub leading corrections and the 
dots stand for terms which vanish faster as Mq — > 0. From Eq. ( |3^ one immediately 
convinces oneself that the anisotropic effects in G{x) vanish for Mq — > as Mq where p is 
a universal critical exponent given by 
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p = 2 + 7]-r]i. 



(35) 



We must notice that the exponent p is not related to the exponent u which characterizes the 
critical behaviour of the "subleading" terms which vanish as Mq, as they are connected to 
different (rotationally-invariant) irrelevant operators. Finally notice that the leading term 
breaking rotational invariance is universal apart from a multiplicative constant, the factor 

Let us now consider the small-momentum limit in which y ^ keeping Mq fixed. In this 
case one can write for / = 0, 2 (or in the special case we have discussed above for / = 0, 1) 

00 

g2i{y,MG)= Y.U2i,m{MG)y^. (36) 

m=0 



By comparing this expansion with Eq. (|30| ) and using Eq. (^), one recognizes that 

Z2i{Mg) = U2l,oiMG), (37) 

and 

C2i,m = hm 7TTT- (38) 

^^G-O U21,o{Mg) 

In the following sections we will use this formula to derive estimates for C2i^m- Indeed 
the functions tt2i,m(^G) can be determined by computing dimensionless invariant ratios of 
moments of G{x; (3): 

q2iAP) = Y.i^'rQ^i{x)G{x; /3). (39) 



It is interesting to notice that the expansion (^^) implies some universality properties 
for some ratios of q2i,m- It is easy to verify that 

p .r,^ _ g0,n(/^)94,m(/^) /^^v 

qO,m{P)q4,n{P) 

is universal for T T^, indeed the constant Z4/Z0 drops out in the ratio. Notice that this 
means that not only Ri^rn,n does not depend on the particular Hamiltonian, but also that 
it is independent of the lattice structure as long as 04 o(x) is the leading operator breaking 
rotational invariance. 



B. Other regular lattices 

All the considerations of the previous subsection can be extended without changes to 
other lattices with cubic symmetry, such as the b.c.c. and the f.c.c. lattices. For other 
Bravais lattices the same general formulae hold, but different multipole combinations will 
appear in the expansion, according to the symmetry of the lattice. In general a larger number 
of multipole combinations with given spin appears when considering lattices with a lower 
symmetry. It is important to notice that in order to have a rotationally-invariant critical 
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limit no multipole Qi{k) with / = 2 should appear in the expansion of the Hamiltonian. 
Thus our considerations apply only to highly symmetric lattices with a tetrahedral or larger 
discrete rotational group. 

As an example of a non-cubic-like lattice let us consider the two-dimensional triangular 
lattice. It is invariant under rotation of tt/3. The relevant multipoles are 

where we have set k^ = \k \ cos 6, ky = \k \ sin 6 and we have assumed one of the generators 
of the lattice to be parallel to the Thus in this case we write 

oo 

k\k)=J2ni{k)eeiiP), (42) 

(=0 

and a similar expression for the expansion of the two-point function. For the triangular 
lattice the first operator which breaks rotational invariance has dimension d + 4. This is a 
consequence of the fact that the triangular lattice has a larger symmetry group with respect 
to the square lattice. We define moments corresponding to TQi{k) by 

hi,m{P) = T.i^Yni{x)G{x; P). (43) 

X 

The arguments given in the previous subsection can be generalized to the triangular lattice 



in a straightforward way. One derives an expansion of the form (34) with p = A + r] — rje, 
T(i{k/MG) and ga^y) substituting QAik/Mc) and g^iy). 



C. Non-Bravais lattices 

Up to now we have considered regular (Bravais) lattices. However other important lattice 
structures are represented by lattices with basis. Particular examples are the honeycomb 
lattice in two dimensions and the diamond lattice in three dimensions. These lattices are 
generically defined by the set of points x such that 



X 



x' + pffp, x' = ^lif]i, (44) 



where f]p is the so-called basis vector joining the two points of the basis, and f/j are the 
generators of the underlying regular lattice. Here p = 0, 1 and U E Z. For the honeycomb 
lattice fji are the generators of a triangular lattice while for the diamond lattice r/j are the 
generators of a f.c.c. lattice. Due to the breaking of translational invariance one distinguishes 
between correlations between points with the same value of p (i.e. points belonging to the 
same regular lattice) and points with different p. In general the components Gppi of the 
two-point correlation function can be written in the form 

G„o(x-y) = G„(.-,) = /|e«'->^^J^ (45) 
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and 



G„.(.-.)^G„(.-.)^/^e.'<^-.'«, (46) 

where the integrals are performed over the Brillouin zone of the corresponding underlying 
regular lattice, Vb being its volume. Gii{x) and therefore A{k,MG) have the symmetries 
of the underlying regular lattice and thus can be expanded as in the first subsection. On 
the other hand, H{k, Mg) does not have the symmetry of the regular lattice, but only 
the reduced symmetry of the full lattice. For the Gaussian model with nearest-neighbor 
interactions defined on the honeycomb and diamond lattices (and also on their (i-dimensional 
generalization), it is easy to realize that, when Mq —>■ 0, 

A{k,MG) — >d[l-\H{k,0)\^]+M^, (47) 

and A(A;,Mg) turns out to be the inverse propagator for the Gaussian theory defined on 
the corresponding regular lattice. In App. |Atj| we present a more detailed analysis of the 
Gaussian theory with nearest-neighbor interactions on the diamond lattice. 

Because of the reduced symmetry, additional multipoles which are not parity-invariant 
appear in the expansion of Hik^Mc)- In the case of the honeycomb lattice the symmetry 
of the triangular lattice is reduced to ^ + Assuming that one of the links leaving a 
site is parallel to the x-axis, one can write 

oo 

H{k,MG) = Y.Tzi{k)h,i{y,MG), (48) 



where we have extended the definition ( |41| ) to include odd multipoles: 



3Z/2 / o, \ 

T,,{k) = {-kY" cos(3/0) = E ( 2 J kf-i^k^f''"^- (49) 



m=0 



The factor i in this equation insures that the functions h^ily, Mg) are real for all /. 
For the diamond lattice one can write 

oo Pi 

H{k,MG)=Y.Y.Qt\k)ht\y,MG). (50) 

/=Op=l 

where Q'f\k) are multipoles constructed from as in the case of the cubic lattices. 

The only difference is that now odd-spin operators are allowed, belonging to the class 

Q2i+3{k) = ikik2k3Q2i{k), (51) 



where we have assumed the natural orientation of the underlying f.c.c. lattice (see App. 

For these lattices, it is not straightforward to make contact with the field-theoretical 
approach. The problem is writing down operators in the effective Hamiltonian that break 
the parity symmetry. These operators must have an odd number of derivatives, but, if they 
are bilinear in a real field 0, they give after integration only boundary terms. The solution 
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to this apparent puzzle comes from the fact that the effective Hamiltonian for models on 
lattices with basis is naturally written down in terms of two fields, defined on the two regular 
sublattices (for the diamond lattice see App. |A3| ). 

As in the regular lattice case, we can associate to the breaking of the parity symmetry a 
universal exponent pp. In principle it can be derived from the critical dimension of the lower- 
dimensional operator breaking this symmetry. From a practical point of view it is simpler 
to consider moments of G{x). For the diamond lattice one defines pp from the behavior, for 
Mq — >• 0, of the odd moments qs^miP), i-e. 

1 ^-3-2m+pp 



Ma""""'''. (52) 



The same formula applies to the honeycomb lattice with the obvious substitutions, go,o 

III. CRITICAL BEHAVIOR OF G{x) AT LOW MOMENTUM 

A. Parametrization of the spherical limit of G{x) at low momentum 

According to the discussion presented in the previous section, in the critical limit multi- 
pole contributions are depressed by powers of Mq, hence for (3 ^ [3c 

^^My). (53) 

where, again, y = k"^ /Mq. As already stated by Eq. (|30|), ^o(y) can be expanded in powers 
of y around ?/ = 0: 

oo 

gQ{y) = l + y + Y.c^y\ (54) 

i=2 



where Cj = Co,i. For generalized Gaussian theories Cj = 0. As discussed in Sec. [11 A| the 



coefficients q of the low-momentum expansion of ^0(2/) can be related to the critical limit 
of appropriate dimensionless ratios of spherical moments 

= qo,j = \x\'^^G{x), (55) 

X 

or of the corresponding weighted moments 

= — . (56) 

It is easy to compute the behavior of fn2j for Mq — > 0: 

m2 ^ 2dMf, (57) 

mi^8d{d + 2){l- C2)Mq\ 

me ^ 4:8d{d + 2){d + 4) (1 - 2c2 + C3) M^^, 
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etc., where d is the lattice dimension. Then, introducing the quantities 



one may compute iii = "Uo^j/mo.o from the following combinations of V2 



■J 



U2 = I- V4, 

U3 = l-2v^ + ve, (59) 

etc... By definition, see Eqs. ( |5B| ) and (0), in the critical limit Ui Ci. 

Another important quantity which characterizes the low- momentum behavior of go{y) is 
the critical limit of the ratio M'^/Mq, 

where M is the mass-gap of the theory, that is the mass determining the long-distance 
exponential behavior of G{x). The value of Sm is related to the negative zero yo of go{y) 
which is closest to the origin by 

yo = -Sm- (61) 

The constant Sm is one in Gaussian models (i.e., when go{y) = 1 + y), like the large- iV limit 
of O(A^) models. 

Let us now consider the relation between the renormalization constant Zg defined at 
zero momentum, 

Zg = xMl = Z.'M'a, (62) 



where Zq has been introduced in Eq. (Pq), and the on-shell renormalization constant Z, 
which is defined by 

OW - (63) 

when k — ^ iM. The mass gap M and the constant Z determine the large-distance behavior 
of G{x); indeed for |x| 00 



M 



d-l 
2 



The critical limit Sz of the ratio Zg/Z is a universal quantity given by 

Z d 

^' = H.-^ = &y^'^y^'^'=''- ^^^^ 
In a Gaussian theory Zq = Z. 
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B. l/AT-expansion 



In the large- limit the difference 

^o(y)-(l + y) (66) 

is depressed by a factor 1/N . It can be derived from the 1/N expansion of the self-energy 
in the continuum formulation. One finds PJT^ 

^o(y) = l + 2/+^0i(2/) + o(^). (67) 

where, for = 3, 
2 

My) = - dz 



^ "'0 arctan 



(Iv^) 



1 y + z + 2^+l \ 1_ ^ y{3 - z) 



4a/P \y + z-2Jyz+l z+l (-2 + 1) 



(68) 



A general discussion of the 0{1/N) correction to go{y) in rf-dimension is presented in App. 0. 
In particular Eq. ( |68D can be derived from Eqs. (|Bl|), ( |B2[ ) and ( [B12| ). The coefficients q of 
the low-momentum expansion of go{y) turn out to be very small. Writing them as 

c'" / 1 \ 



one obtains 



= -0.00444860..., (70) 
4^) = 0.000134410..., 
= -0.0000065805..., 
= 0.0000004003..., 

etc.. For sufficiently large A^ we then expect 

Ci < C2 < 1 for i> 3. (71) 

As we shall see from the analysis of the strong-coupling expansion of G{x), the pattern (|7lD 
is verified also at low values of A^. 

The ratio 5*^ = M^/Mq is obtained by evaluating the negative zero yo of go{y) closest 
to the origin: 

Sm = -yo = l + ^M-l) + O (^) , (72) 

where 

0i(-l) = -0.00459002.... (73) 
Moreover using Eq. (p3|), one finds 
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S, = l + l^;(-l) + 0(ij) (74) 

where 

= 0.00932894.... (75) 

As expected from the relations ( |7TD among the coefficients Cj, a comparison with Eqs. ( |70D 
shows that the non-Gaussian corrections to Sm and Sz are essentially determined by the 
term proportional to (A;^)^ in G~^{k), through the approximate relations 

Sm^I + C2, (76) 
5z^l-2c2, (77) 

with corrections of 0(03). 



C. g-expansion in three dimensions 

Another approach to the study of the critical behavior in the symmetric phase of 0{N) 
models is based on the so-called (^-expansion, the perturbative expansion at fixed dimension 
d = 3 for the corresponding 0^ continuum formulation |T^. The perturbative series which 



are obtained in this way are asymptotic; nonetheless accurate results can be obtained using 
a Borel transformation and a conformal mapping which take into account their large-order 



behavior. As general references on this method see for instance Refs. and [Q. This 
technique has led to very precise estimates of the critical exponents. 

Starting from the continuum action (P), one renormalizes the theory at zero momen- 
tum using the following renormalization conditions for the irreducible two- and four-point 
correlation functions of the field 0: 

r(^Hp)a, = ^c'rg)(p)5.,, (78) 

r(^)(0, 0,0,0)^^^5 = -Za^^McS^p^s, (79) 

where 

T^r\p) = MI, + p' + 0{p'), (80) 

and Sai3^s = Sa/sS-ys + SaySps + SasSj3y. When Mg — > the renormalized coupling constant is 
driven toward an infrared stable zero g* of the /5-function 

The universal function go{y) is related to the renormalized function 

f{g,y)^MfT^i\k) (82) 

by 
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go{y) = lim 



(83) 



We computed the first three non-trivial orders of the non-Gaussian corrections to go{y)- A 
calculation up to four loops gave 



f{g,y) = l + y + fZ^Z^ 



-4^4^2 N + 2 



{N + 



2v2^ N + 2 , , , -3v3v3/2 N + 2 

9^^jNTW^^iy) + 9 z,z^ ^^^^ 

{N + 2) , , {N^ + 6N + 20) , , (5iV + 22) , / 



where g is the rescaled coupling 



{N + 



(N + t 
487r 



(84) 

+ Oif), 



Zg is the renormalization constant of the coupling (defined by go = McgZg) 

2(41A^+ 190)' 



1+^ + 



1 - 



27(A^ + 8)2 

and Zg is the zero-momentum renormalization of the field 



r + o{g 



27(A^ + 8) 



(85) 



(86) 



(87) 



A simple derivation of the two and three- loop functions ^2{y) and ^z{y) is presented in 
App. (cfr. Eqs. ( P19|) ). In particular using the results of Refs. [|l5],|l6l one finds 



1 X arctan(./w/3) 4 
if2{y) = 41n(l + ly) + 24 " ^ " 



(88) 



We shall not report the expressions of the four-loop functions fijiy) because they are not 
very illuminating. 

The coefficients of the low-momentum expansion can be easily obtained from Eq. (|85D 
by calculating the zero-momentum derivatives of the functions (pn,j{y)- We write 



iN + 



{N + 



2h? + m f 



N + 2 

{N + sy 



2h 



(2) 



8(7iV + 32) 
3(iV + 8)2 



+?>hf^ + h 



(4,1) {N + 2) ^ ^(,^,) (iV2 + 6Ar + 20) , ^(4,3) (SAT + 22) 1^ 



(AT + 8)2 ' (AT + 8)2 

where we have introduced the coefficients 



hT 



>t + 0{f), (89) 



(A^ + 



h 



1 d' 
i\ dy 



J^n,j{y)\y=0- 



(90) 



In Table | we report the numerical values of h^^'^'^ for i < 5. 

By evaluating the zero of goiy) closest to the origin, one obtains 
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s 



M 



(AT + 8)2 
which numerically leads to 

2 N + 2 



Sm - 1 



-9 



(AT + 8)' 



-0.00521783 



(91) 



1 + ^ X 0.054182 + 0(^^ 



Moreover 



(AT + 8 

which numerically leads to 



A^ + 2 
(AT + 8)' 



r0.0107349 



1+^ X 0.041829 + 0(^2 



(92) 



(93) 



(94) 



A comparison of the (^-expansions of q, 5*^/ and Sz shows that the approximate relations 
([TBI) and ( |77D are valid for all values of A^ and not only for A^ — > oo as shown in the previous 
subsection. 

In order to get quantitative estimates, one must perform a resummation of the series 
and then evaluate it at the fixed-point value of the coupling g* . Although the terms of the 
(yf-expansion we have calculated are only three for q and two for Sm and Sz-, we have tried 
to extract quantitative estimates taking into account also the following facts: 

(i) The (yf-expansion is Borel-summable [T^ (see also e.g. Refs. p| and [|14| for a discussion 
of this issue), and the singularity closest to the origin of the Borel transform (corresponding 
to the rescaled coupling g) is known |18[: hg = —0.75189774 x (A^ + 8). 

(ii) The fixed point value g* of g has been accurately determined by analyzing a much 
longer expansion (to 0{g'^)) of the corresponding /5-function ||T9|^2|. Indipendent and con- 
sistent estimates of g* have been obtained by other approaches, such as strong-coupling 
expansion of lattice non-linear O(A^) a models [^,2^ (for A^ = 1 see also Refs. p5-p8[]), 
and Monte Carlo lattice simulations (only data for A^ = 1 are available p9|-p^). 

We have followed the procedure described in Ref. |3^ (see also Ref. 0), where the per- 
turbative expansion in powers of g is summed using a Borel transformation and a conformal 
mapping which takes into account its large-order behavior. We transform the series 



R{9) 



k=0 



into 



R{g) = Y.Bk e"'u{gtfdt 



(95) 



(96) 



where 



(97) 
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The coefficients are determined by comparing the expansion in powers of g of the r.h.s. of 
Eq. (^) with the original expansion. Since the (^-series of Cj, S'a^ — 1 and Sz — l, have the form 
^{9) = 9'^J2k=o(^i9\ may apply the resummation method to either R{g) or R{g)/g'^. 
In Table we present results for both choices. Following the suggestions of Ref. |33| we 
also tried more refined resummations, changing formula (^) to weaken the singularity of 
the Borel transform. We did not find any significant difference. In our calculations we used 
the estimates of g* obtained from the analysis of the /5-function by JT^pn| , P^ . They are 
reported in Table 0. For small values of slightly lower values of g* were computed in Ref. 
taking into account the non-analiticity of the /3-function at the critical point 



This 

difference is however too small to be quantitatively relevant in our calculations. 

It is difficult to estimate the uncertainty of the results: the fluctuations of the results 
with respect to the method we used to resum the perturbative series indicate an error of 
^ 20% on Ci and Sm for small values of N. As increases the estimates become more 
stable. The final results are in good agreement with the estimates by other methods. 



D. e-expansion 

The universal function go{y) can be computed perturbatively in e = 4 — d using the 
continuum 0^ theory The leading order is simply go{y) = 1 + y. The ffist correction 
appears at order and was computed by Fisher and Aharony ||^. We have extended the 
series, calculating the O(e^) term, obtaining 



goiy) = l + y + e' 



N + 2 



1 + 



6(3A^ + 14) 
(AT + 8)2 



N + 2 

(iv + sy 



.i;,(y) + Oie'), 



(98) 



where 



My) = 2 ^z{l + \z)\n (^^l + \z + Iv^) h{y, z), (99) 
h{y, z) = + jr-^ + 2^(^+^ + '^~ ,Jl + 2y + 2z + y^- 2yz + z^^ . 



l + z 



l + z) 



We do not report the explicit expression of ipsiy) because it is not very illuminating. It can 
however be obtained from Eqs. ( |B15| ), ( |B16| ) and ( |B18| ) of App. ^, where we show how to 
derive the functions ijj2{y) and ip3{y) from the 0(1/A^) calculation of go{y) in d dimensions. 
Setting 



A^ + 2 



one finds 



C2 = -0.00752024... x 
C3 = 0.000191931... X 



1 + e 




0.249301... +0(e^ 



(100) 



(101) 



1 + e 



/6(3Ar+ 14) 
1^ (AT + 8)2 



0.130607... 
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C4 = -0.0000081420... x 



l + ,f5M±ll)_0.003053...Vo(a 



C5 = 0.0000004391... x 



1 + e 



'6(3A^+ 14) 



+ 0.117278... +0(6^ 



etc.... 

One also obtains 



S 



M 



1 + 



N + 2 

(ivTsy 



6(3Ar + 14) ■ 
^ (iV + 8)2 ' 



M-l) + Oie'), (102) 



where 



^2(-l) = -0.00772078..., tpsi-l) = 0.00189984. 



and 



where 



, N + 2 

{N + sy 



6(3Ar+ 14) ■ 
^ (iV + 8)2 ' 



^^(-l)+6^ 



N + 2 



{N + 



^^(-1) = 0.0156512. 



^P'^{-l) = -0.0038246.... 



(103) 



(104) 



(105) 



In order to get quantitative estimates from the perturbative e-expansion, one should first 
resum the series and then evaluate the resulting expression at e = 1. Usually resumma- 
tions are performed assuming the Borel summability of the e-series. As in the case of the 
^f-expansion, a considerable improvement is obtained if one uses the knowledge of the singu- 
larity of the Borel transform [|18[, bg = —{N + 8)/3. We have used the procedure described 
in the previous subsection. Again, since the e-series of Cj, 5*^/ — 1 and Sz — I, have the 
form R{e) = =o^j^*5 we applied the resummation method to i?(e) and to R{e)/e^. In 
Table |^ we present results for both choices. Since we use a series with only two terms the 
estimates are not very precise as the large difference between the results obtained with the 
two methods indicates. 

One can also try to get estimates for two-dimensional 0{N) models, i.e. for e = 2. 
By resumming the series of C2(e) and S'A/(e), we find: C2 = —0.0010 and Sm = 0.9989 for 
N = 1, which must be compared with the exact results C2 = —0.000793... and Sm = 



0.999196...; Ca = -0.0013 and Sm = 0.9987 for = 3, to be compared with the strong- 
coupling results C2 = —0.0012(2) and Sm = 0.9987(2). In both cases the agreement 
is satisfactory. Instead, when resumming the series divided by the agreement is poorer. 
We find C2 = -0.0026 and Sm = 0.9973 for A^ = 1 and C2 = -0.0028 and Sm = 0.9971 
for A^ = 3. A posteriori, it thus appears that the estimates obtained from the resummation 
of the complete series R{e) are more reliable. This is confirmed by the three-dimensional 
analysis where the estimates obtained by considering -R(e) are those which are in better 
agreement with the strong-coupling and (yf-expansion estimates. 

For quantities which are exactly known in two dimensions, one can modify the e-series 
to obtain a new expansion which gives the exact value for e = 2. This can achieved |^ by 
defining for a generic observable R, with e-series -R(e), 
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R{e) 



R{e) - R{e = 2) 



(106) 



and then rewriting 



R{e) = i?^^"'^* + (2 - e)R{e) 



(107) 



where R'^^°-'^^ is the exactly known value of R in two dimensions. In other words one subtracts 
to the original series its value for e = 2 and then sums the exact two-dimensional result. 
Estimates for e = 1 are obtained by resumming the perturbative expansion of -R(e) and 
computing -R(l). We applied the method to the Ising model, i.e. = 1, since in this case, 
the coefficients Cj and Sm are exactly known [^. As before, one can also apply the same 
procedure to R{e)/e^, defining 



Me) 



R{e) R{t = 2) 



(108) 



and then writing 



R^ie) 



-R 



exact 



+ e\2-e)R2ie). 



(109) 



From a conceptual point of view, 
(|107|) gives the exact value for d = 



Eq. (|109D appears preferable to Eq. (|107|) . Indeed Eq. 
: 2, but it does not reproduce the correct value for d = 4. 
Eq. (|109|) instead predicts correctly -R2(e) ~ O(e^) for e — > 0. In any case we report 
the results obtained with both methods in Table [V|. They are referred to as "improved" -e 
expansion. The estimates are in good agreement with the other results. Notice also that the 
large discrepancy between the two different resummations of the unconstrained e-expansion 
is here significantly reduced. 



E. Strong-coupling analysis 

In this subsection we evaluate some of the quantities introduced in Sec. |III A| by analyzing 
the strong-coupling expansion of the two-point function G{x) in lattice O(A^) non-linear a 
models. We recall that non-linear a models and 0'^ models possessing the same internal 
symmetry 0{N) describe the same critical behavior, thus giving rise to the same universal 
two-point function go{y) in the critical limit Mq 0. 

By employing a character-like approach [^, we have calculated the strong-coupling 
expansion of G{x) up to 15th order on the cubic lattice and 21st order on the diamond 
lattice for the corresponding nearest-neighbor formulations. In App. we present the 15th 
order strong-coupling expansion of G{x) on the cubic lattice. In App. ^we report the 21st 
order strong-coupling series of the magnetic susceptibility and of the second moment oiG{x) 
on the diamond lattice for = 1, 2, 3. 

We mention that longer strong-coupling series, up to 21st order, of the lowest moments 
of G{x) on the cubic and b.c.c. lattices have been recently calculated by a linked-cluster 
expansion technique, and an updated analysis of the critical exponents 7 and u has been 
presented PD| . For A^ = even longer series have been calculated for various lattices pil-HH . 
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In our strong-coupling analysis, we took special care in devising improved estimators 
for the physical quantities Cj and Sm, because better estimators can greatly improve the 
stability of the extrapolation to the critical point. Our search for optimal estimators was 
guided by the large- limit of lattice O(A^) a models. 

In the large- limit of O(A^) a models on the cubic lattice the following exact relations 
hold in the high-temperature phase, i.e. for (3 < jSc, 

uriMc) ^u, = ^M^, (110) 

etc., which vanish for T Tc, i.e. for Mq — > 0, leading to the expected result q = 0. 
Similarly on the diamond lattice one obtains 



1 1 _|_ 3/M'2 



etc... The above formulae can be derived from the large- limit of the two-point function 
on the cubic and diamond lattice given in App. ^ 
We introduce the following quantities 

u, = u,-u^{Mg), (112) 

whose limits for T ^ are still c,. At N = oo Ui are optimal estimators of q, since 

u,{p) =u* = Ci = (113) 

for P < Pc, i-e. off-critical corrections are absent. It turns out that the use of Ui, beside 
improving the estimates for large values of A^, leads also to more precise estimates of Cj at 
low values of A^. Strong-coupling series of Ui can be easily obtained from the strong-coupling 
expansion of G{x). 

On the lattice, in the absence of a strict rotational invariance, one may actually define 
different estimators of the mass-gap having the same critical limit. On the cubic lattice 
one may consider fi obtained by the long-distance behavior of the side wall-wall correlation 
constructed with G{x), or equivalently the solution of the equation 

G"^(i/i,0,0) = 0. (114) 

In view of a strong-coupling analysis, it is convenient to use another estimator of the mass- 
gap derived from /i p4|,|6|: 

= 2 (cosh/i - 1) , (115) 

which has an ordinary strong-coupling expansion (/i has a singular strong-coupling expan- 
sion, starting with — ln/3). One can easily check that M^/n — 1 in the critical limit. A 
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similar quantity M| can be defined on the diamond lattice, as shown in App. |A 3| (cfr. 
Eq. ( A25|) ). One may then consider the dimensionless ratios M'^/Mq and M^/Mq respec- 



tively on the cubic and diamond lattices, and evaluate their fixed-point limit Sm-, which by 
universality must be the same. 

In order to determine the coefficients C2 and C3 of the low- momentum expansion of go{y) 
and the mass-ratio 5*^/, we analyzed the strong-coupling series of U2 and (defined in 
Eq. ( p.l2| )), and those of the ratios M'^/Mq and M^/Mq respectively on the cubic and 
diamond lattice. 

On the cubic lattice the available series of U2, and M^/M^ — 1 are respectively of 
the form 13^ Yh=q (^il3\ J2t=o o-iP', and SLo ^j/^*! except for = 1 where they are of 
the form (3^ Y.I=q ai(3\ (3^ T.I=q ai(3\ and (3^ T.I=q ai(5\ These series can be derived from the 
strong-coupling expansion of G{x) presented in App. 0. On the diamond lattice the available 
series of ^2, and M^/Mq — 1 are respectively of the form J^l^o J2ito and 

P^ELo O'iP''^ except for A^ = 1 where they are of the form P^ElioaiP\ P'^ EHo (^iP' ^ and 
P^Elo^iP' 

We constructed various types of approximants to the above series, such as Pade approx- 
imants (PA's), Dlog-Pade approximants (DPA's) and first-order inhomogeneous integral 
approximants (lA's) We then evaluated them at the critical point jSc in order to obtain 
an estimate of the corresponding fixed point value. For the cubic lattice and most values of 
A^, jSc is available in the literature from strong-coupling and numerical Monte Carlo studies 
(see for example Refs. p3|j40|j4^j46| -^ ) . When /3c was not known (as in the case of diamond 



lattice models for A^ > 0), we estimated it by performing a lA analysis of the strong-coupling 
series of the magnetic susceptibility x = J2xG{x). In our analysis errors due to the uncer- 
tainty on the value of Pc turned out negligible. The values of Pc used in our calculations are 
reported in Table fT|. 

In the analysis of a series such as A = P"^ J27=o (^iP\ 'we constructed approximants to 
the nth order series P~''^A = Z^iLo'^*/^*' then derived the original quantity from them. 
Given a nth order series, we considered the following quasi-diagonal approximants: 

(i) [l/m] PA's with 



(ii) [l/m] DPA's with 



I + m > n — 2, 

Ti 

l,m>--2; (116) 



l + m + l>n-2, 

71 

l,m>--2; (117) 

(iii) [m/l/k] lA's with 

■m + l + k + 2 = n, 

[{n - 2)/3j -l<m,l,k< \{n- 2)/3 \ + 1. (118) 

In the case (i) and (ii), /, m are the orders of the polynomials respectively in the numerator 
and denominator of the PA of the series at hand, or of the series of the logarithmic derivative 
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in the case of DPA. In the case of integral approximants, m,l,k are the orders of the 
polynomial Qm, Pi and Rk defined by the first-order linear differential equation 

Qm{x)f{x) + Pi{x)f{x) + R,{x) = O , (119) 

whose solution provides an approximant of the series at hand. 

In Table |ITT| we show some details of the above-described analysis for the series on the 
cubic lattice and for selected values of A^, including those physically interesting. We report 
there various estimates as obtained from the corresponding plain series, and the resumma- 
tions by PA's, DPA's and lA's. On the cubic lattice, for A^ < 8, the plain series of U2 already 
provides a good estimate of C2, indeed the values at Pc of the series using i terms appear 
to oscillate around an approximately constant value. In these cases as an estimate of C2 we 
can take the average of the values of U2 at f3c obtained using n and n — 1 terms (the errors 
reported in Table |IT| are related to their difference). For the C2 corresponding to other values 
of A^ and for C3 and Sm, such oscillations are not observed but the results from the plain 
series are close to their resummations. In these cases in Table PTT| we report just the value at 

using all available terms of the series. From the PA, DPA, lA analysis we take the average 
of the values at /3c of the non-defective approximants using all available terms of the series. 
As an indicative error from each analysis we consider the square root of the variance around 
the estimate of the results from all non-defective approximants specified above. PA's, DPA's 
and lA's are considered defective when they present spurious singularities close to the real 
axis for Re/5 ^ Pc- More precisely we considered defective the approximants with spurious 
singularities located in the region < Re/5 < 1.1 x /5c (sometimes < Re/? < 1.2 x jSc) and 
|Im/5| < 0.5 X /3c- Anyway, our final results turned out to be quite insensitive to this choice. 
Most of the approximants we considered turned out non-defective. Similarly in Table |IV| we 
report results of the analysis of the series on the diamond lattice. In this case we do not 
report estimates from the plain series because they appear not to be reliable and very far 
from their resummations. 

Table summarizes our calculations. The final estimates of C2, C3 and Sm reported in 
Table synthetize the results from all the analyses we performed, and the reported errors 
are a rough estimate of the uncertainty. Final results are rather accurate taking into account 
the smallness of the effect we are looking at. Universality among the cubic and diamond 
lattices is in all cases well verified and gives further support to our final estimates. Our 
results are in good agreement with the estimates obtained from the g- and 1/A^ expansion. 
Only for C2 and small values of A^ one notices a small discrepancy, probably due to confluent 
singularities, which are not properly handled by the approximants we considered. We also 



tried different resummation methods |50| which take into account the subleading corrections. 
However, in this case, most of the approximants turned out to be defective and no reliable 
estimate could be obtained. 

Our strong-coupling analysis represents a substantial improvement with respect to earlier 
results reported in Ref. ||^ for the Ising model, and obtained by an analysis of the strong- 
coupling series calculated in Refs. fllSTI: C2 = —5.5(1.5) x 10~^, C3 = 0.05(2) x 10""^ on 
the cubic lattice, and C2 = —7.1(1.5) x lO"'' and C3 = 0.09(3) x 10"^ on the b.c.c. lattice. 



Other strong-coupling results can be found in Ref. |Q . Our analysis achieves a considerable 



improvement with respect to such earlier works essentially for two reasons: we use longer 
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series and improved estimators, see Eq. ( |112|) , which allow a more stable extrapolation to 
the critical limit. Estimates from the analysis of the strong-coupling series of the standard 
variables Ui, defined in Eq. (^9]), are much less precise, although consistent with those 
obtained from Uj. 



F. Summary 

We have studied the low-momentum behavior of the two-point function in the criti- 
cal limit by considering several approaches: 1/A^-expansion, (^-expansion, e-expansion and 
strong-coupling expansion. A summary of our results can be found in Table 0. 

From the analysis of our strong-coupling series we have obtained quite accurate estimates 
of the coefficients C2, C3 of the low momentum expansion (|5^). Asymptotic large- formulae 
and d?^ ) are clearly approached by our strong-coupling results, but only at rather 
large values of A^. The same behavior was already observed for other quantities like critical 
exponents and the zero-momentum renormalized four-point coupling [^. We have also 
computed the universal function go{y) in the ^f-expansion in fixed dimension to order 0{g^) 
and in e-expansion to order 0{e^). The corresponding estimates of C2, C3 and Sm are in 
reasonable agreement with the strong-coupling results. 

For all values of A^ the coefficients C2 and C3 turn out very small and the pattern ([71|) is 
verified. Furthermore relation ([75]) is satisfied within the precision of our analysis. 

The few existing Monte Carlo results for the low-momentum behavior of the two-point 
Green's function are consistent with our determinations but are by far less precise. Using 
Refs. [p2|-[5^ one estimates C2 = —13(17) x 10~^ for self- avoiding walks, which correspond 



to A^ = 0. In Ref. I^Sf the authors give a bound on VSm for the Ising model (A^ = 1), from 
which —1.2 X 10"'^ < Sm — 1 < 0, which must be compared with our estimate Sm — 1 = 
—2.5(5) X 10""^. Monte Carlo simulations of the XY model (A^ = 2) shows that Sm — 1 
within 0.1% which is consistent with our strong-coupling result ^a/ — 1 — —3.5(5) x 10""^. 

We can conclude that in the critical region of the symmetric phase the two-point Green's 
function for all A^ from zero to infinity is almost Gaussian in a large region around fc^ = 0, 
i.e., \k'^ /Mq\ % 1. The small corrections to Gaussian behavior are dominated by the {k'^Y 
term in the expansion of the inverse propagator. Via the relation (|I|) such low-momentum 
behavior could be probed by scattering esperiments by observing the low-angle variation of 
intensity. A similar low-momentum behavior of the two-point correlation function has been 
found in two-dimensional 0{N) models [p6| , p7| , p6| . Substantial differences from Gaussian 



behavior appear at sufficiently large momenta, where G{k) behaves as 1/fc^ ^ with f] ^ 
(although rj is rather small: rj ~ 0.03 for < A^ < 3). 



IV. ANISOTROPY OF G{x) AT LOW MOMENTUM AND IN THE CRITICAL 

REGIME 

In this Section we will study anisotropic effects on the two-point function due to the 
lattice structure. We will mainly consider three-dimensional lattices with cubic symmetry. 
However, whenever possible, we will give expressions for general d-dimensional lattices with 
hypercubic symmetry, so that one can recover the results for the square lattice and compare 
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with perturbative series in d = A — e. We will also comment briefly and present some results 
for the triangular, honeycomb and diamond lattices. 



A. Notations 



In the following subsections we will compute the exponent p = 2 + rj — rj^ deflned in Eq. 
(^). It can be derived directly from Eq. ([3l| ) or Eq. (|3^) or by studying the weighted 
moments q^^j = q4,j/mo where is deflned in Eq. (|^) and mo = X- Indeed for Mq — > 0, 



q,^^ ~ M^'-'™+^ (120) 

We will also compute the universal function gi{y). In particular we will be interested in the 
flrst terms of its expansion in powers of y around y = 0: 

g,{y) = l + Y.d^y\ (121) 

where di = C4^i (cfr. Eq. (|30D). The coefficients di can be easily obtained from the expressions 
of the moments q^^m- For Mq 0, we find 



4(d + 8) (l-^dAuf, 



94,1 
^4,0 

^ — > 2A{d + 8)(rf + 10) fl - Ui - |c2 + \d2) (122) 

?4,0 ^ ^ 



and so on. From (|122|) it is easy to derive expressions for = M4,i/M4,o whose critical limit 
is di. In particular 

Ml qAi , , 

n = 2 - , , ^ , 123 



B. Breaking of rotational invariance in the large- AT limit 

In the large- limit lattice O(A^) models become massive Gaussian theories which can 
be solved exactly. If one considers theories deflned on Bravais lattices one has in the large- iV 
limit 

G"\k) =cp{k^ + Ml), (124) 

where k is deflned by Eq. (^. The relation between Mq and (3 is given by the gap 
equation. The constant c is lattice-dependent and will not play any role in the discussion. 
Speciflc examples are the nearest-neighbor Hamiltonians of the form (^ for which we have 
collected in Appendix ^ general expressions for various three-dimensional lattices, the cubic, 
diamond and f.c.c. lattice. Analogous formulae for some two-dimensional lattices, the 

1-2 



square, triangular, and honeycomb lattice can be found in Ref. ||5H]. The function k has the 
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properties mentioned at the beginning of Sect. P A| and a multipole expansion of the type 



(10) for lattices with cubic symmetry. For other Bravais lattices the only difference is the 
presence of different multipole combinations. Considering first lattices with (hyper)-cubic 
symmetry, from Eqs. ([I0| ) and (|18D, we find for Mq — >■ 



G-\k) = cl3Ml \l+y+ M2 (62,02/' + e.^oQ^ik/ Mg)) ..]. (125) 



Comparing with Eqs. (^) and ( P^ we get immediately p = 2 and gi{y) = 1, i.e. di = for 
all i ^ 0. 

In the large- iV limit one can easily verify the universality properties of the ratios defined 
in Eq. (|^. Indeed for generic Hamiltonians in the critical limit Mg —>■ (keeping the 
dimension of the lattice d generic) we have 



/m—l 



2'"m! [l[{d + 2i)] Mg^"^, (126) 



(l4,r. 



, i=0 

/m—l 



2'^(m + l)! n(^ + 8 + 2^)p^G^™' (127) 
and 

q,,o ^ -64,0 ' m^\ (128) 

Notice that the only dependence on the specific Hamiltonian is in the expression of q^ Q. 
(Exact expressions for some of these quantities are reported for the theory with nearest- 
neighbor interactions on the cubic, diamond and f.c.c. lattice in Table |V1| and on the square 
lattice in Table [VI H ) . Universality is then a straightforward consequence of the independence 



of the ratio (|127|) from 64,9. It should also be noticed that 

A4,m = ~ Ml (129) 

This shows that, as expected, anisotropic moments are suppressed by two powers of Mg in 
agreement with the prediction p = 2. We stress that the universality of Ri,m,n is due the 
fact that there is only one leading irrelevant operator breaking rotational invariance. 

It is interesting to notice that such a universality does not hold for moments gg ^ (or for 
^2/ ra higher values of /) because of the mixings we have mentioned in Sec. |II A| . For ^ 
we have for T ^ Tc 

^6,- 2™(m + 1)! f 1 + ^7^1 ill id + 12 + 2^)) M^'"^, (130) 



gg ,0 V 66,0 d+12 



=0 



which depends on 65,0 and e^Q, a consequence of the fact that Qii^k)"^ contains a term of the 
form k'^Q^lk). Thus ratios of the form ( ^OD built with are not universal. 

Let us now consider the diamond lattice. In this case not only rotational invariance is 
broken, but also the parity symmetry. As the leading anisotropic operator is 04,0(0;) the 
behaviour of the leading anisotropic corrections is identical to that we have discussed above. 
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Therefore p = 2 also in this case. Moreover the invariant ratios -R4,m,n are identical for the 
diamond lattice and for the other Bravais lattices with cubic symmetry. Eq. ( 127 ) is exact 
for the diamond lattice as well. 

To discuss parity-breaking effects we must consider odd moments of G{x). In particular 
one finds that, for Mq — > 0, 

mo 6\/3 

where 

93,0 = '^(^' (132) 

Thus parity-breaking effects vanish as Mq,, i.e. Pp = 3, faster than the anisotropic effects 
we have considered previously. 

Finally let us consider lattices which do not have cubic invariance, such as the trian- 
gular and the honeycomb one. In Table [VI I| we report the large-A^ limit of some of the 
lowest spherical and non-spherical moments of G{x) for the models with nearest-neighbor 
interactions. 

For the triangular lattice one should consider the multipole expansion (|4^) . In this case 
the leading term breaking rotational invariance is proportional to TQ[k) and thus we have 
p = 4. This is indeed confirmed by the fact that, for Mq —>■ 0, 

^6,™ = ~ M'a, (133) 

^6+2m 



where tg.m = ^6,m/''^o and tg^m is defined in Eq. (|3]). As in the cubic case, it is immediate to 
verify the universality of ratios of the form given in Eq. (^) with tg ,„ instead of q^^m, which 
is a consequence of the uniqueness of the leading operator breaking rotational invariance. 
Universality follows from the fact that, for T — > T^, 



te^m 2^'"(m + l)!(m + 5)! 
~^ 5! 



M^^™, (134) 



independently of the specific Hamiltonian. 

For the honeycomb lattice one must also consider the breaking of parity. Considering 
the odd moment 



(cfr. Eq. (|^)) one finds 



h,o = J2i^'-5y^x)Gix,y), (135) 



ho - ^ — ^. (136) 
mo 2 



Thus, as in the diamond case, parity breaking effects vanish as Mq, i.e. pp = 3. 
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C. Analysis to order \/N 



In the previous subsection we computed the exponent p for oo for lattices with 

cubic symmetry, finding p = 2. Now we want to compute the corrections, i.e. the value 
oi a = = rj — rj4^ (cfr. Eq. (p5|)), which is the anomalous dimension of the operator 04^0(2;)- 
More generally we can compute the exponents ri2i defined in Eq. ( pi]) for arbitrary /. Notice 
that in this way we will also obtain the correction to p for the triangular lattice which 
depends on rj^. 

In d dimensions, we consider the following representation of the inverse two-point function 
where the 0{1/N) correction has been included 



r ^ t + r 'z-Mg + 1 / ^A(p) ==i-^ -^]. (137) 



— 2 

Here k is the inverse lattice propagator defined in Eq. 

^1 



2 dkl 
and 



A;u=0 



(138) 



2 7 (27r)^(^2^M2)(g2 + M2)' 



The following statements can be checked explicitly in Eq. ( |137| ) and hold to all orders 
of the 1/A^ expansion: (i) in the limit Mq —>■ the function G~^{k, Mq) is spherically 
symmetric (i.e. it depends only on y = k'^/MQ, apart from an overall factor); (ii) the only 
non-spherically symmetric contribution that may appear in G~^{k, Mg) to 0{Mq) can be 
reduced to a spherically symmetric function multiplied by Qi{k). These statements are 
simply a consequence of applying the discrete and continuous symmetry properties to all 
integrals appearing in the asymptotic expansion in Mq of the relevant Feynman integrals. 
They prove to all orders in 1/A^ the validity of the expansion (Q). 

To compute the anomalous dimension 772; to order 1/A^ we will use the trick we explained 
in Sec. [11 A] . If one considers a particular Hamiltonian such that g2i{yi Mq) = for < / < I, 



then G~^{k) has an expansion of the form (|3^ with 774 — > r^gl gA{y) —>■ ^2l(^)- 
1/A^-expansion, to order 1/A^ this can be achieved by considering Hamiltonians such that, 
for k 0, (to simplify the notation from now on we write / instead of I) 

■f = rk^i 0{e'^^) (140) 

where fc^' = The limit Mq — ^ can then be easily obtained by evaluating massless 

continuum integrals, and taking the contribution proportional to r, which is the only term 
relevant to our computation. In this limit we obtain 
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where 



{27rY [{q + + r{q + pf^] [q^ + rq 



2n 



d'^q 



— r 



d'^q 



^21 



2 J {2nY q^{q + py 



^ (47r)2r(c;-2) 



{27TY{q^y{p + qy 



(141) 



(142) 



Bi = iA- d) 



r (f + 2/ - 2) r (c/ - 2) 

r(d + 2/-3)r(f-i 



(143) 



and we have discarded rotationally invariant terms proportional to r, since they will not 
contribute to the final result. 

We must now identify the singular contribution in the limiting form of Eq.( |137| ): 



d'^p 



N J {2t^Y 
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d'^p Ao(p) 



[,2/ 1 /■ dS 
T^nJ(2^ 



l-rBi^— 
d'^p Ao(p) 



{p + + r{p + kf^ 

p2« _|_ J^yi 

Bi— - 

pz 



-1 



(p + ky 



{2'KY{p+kY \ NJ {2TTY{p+ky 

^k'(^l- ^r^i In k^ + r^' (^1 - ^7]2i,i In k^ . (144) 
The coefficients ri\ and ri2i,i are related to the 1/A^ expansion of the exponents r) and ri2i- 



and 



mi 



V2l,l 

N 



By simple manipulations one obtains 



4r{d-2) 



r(2-f)r(f-2)r(f-i)r(f + 1)^ 



(145) 



(146) 



(147) 



and 



V2l,l 



d{d-2) 



(rf-2 + 4/)(d-4 + 4/) 



1 + 2 



T{2l + l)T{d-2) 
r(2/ + rf-3) 



^1- 



(148) 
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Notice the following properties: 



V2,i = Vi for d, (149) 

V2i,i Vi for d^2, (150) 
3 

^2^1 = _ for = 3, (151) 
3 

r]2i,i -n—, 7T?7i for d A. (152) 



/(2/ + 1) 



Therefore for c? = 3 we find 



37r 



^1 = (153) 
which agrees with the known result, and 
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^ = ^-^4 = 7TT-^ + — ). (154) 



For ?7 also the 0{1/N'^) and 0(1/A^^) corrections are known |^^. The coefficient cxi of the 
1/A^ expansion is very small. Thus, at least for sufficiently large, say N ^ 8, where the 
1/A^ expansion is known to work reasonably well, corrections to the Gaussian value of p are 
very small. In two dimensions, and to 0{1/N), there are no corrections to the Gaussian value 
because of Eq. ( |150| ), i.e. the first coefficient of the expansion of the anomalous dimension is 



zero to 0{1/N). One might only observe (suppressed) logarithmic corrections to canonical 
scaling for all /. It is easy to check in perturbation theory that this holds exactly for all 
N>3. 

Alternatively the exponents ri2i could have been computed from Eqs. ( jn\ ) and (PT|). 

The computation of the universal function ^4(2/) is particularly involved and is presented 
in App. |C[ In this Section we will only give the values of the coefficients di which appear in 
the low- momentum expansion (cfr. Eq. ( p.21| )). We found 

0.00206468... ^ / 1 \ ,,,,, 



0.00007378... ^ / 1 \ 
0.00000424... ^/ 1 



etc. 



D. gi-expansion analysis 

The critical exponent a and the scaling function g4^{y) can also be evaluated in the g- 
expansion. For this purpose we calculated the one-particle irreducible two-point function 
Toi{k, Mg) defined in Eq. ([27|). By a three-loop calculation one finds for the bare correlation 
function 
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ToAK Mg) = Q,{k) + gl^^^J2{K Mg) 



where 



and 



A{p, m) 



d^p Qiik — p)A{j),m) 

(27r)3 [(k-p)2^m'^f ' 

c^^P Q^ik — p)A{p,my 

(27r)3 p_p)2 + ^2]2 ' 

(i^p m)AQ(p, m) 
(27r)3 (A;-p)2 + m2 ' 



1 



1 p 

— ctrctj cLll 

(27r)3 [g2 -|- 77^2] [(g -|- p)2 _|_ 2m' 



rf3g 



Q4(g) 



(27r)3 [g2^_^2]2[(g_^p)2_^^2]- 

By appropriately renormalizing To4^{k, Mg) at /c = 0, one derives 
N + 2 AOn^ 



(156) 

(157) 
(158) 
(159) 

(160) 
(161) 



-Q^ {d/de)j2{k,l)\k=o 



(162) 



+9 \^r^L^^f Q^ i9/dk) [3J2{k, 1) - 47r J3,i(A;, 1) - 16n Js,2{k, 1)],^^ + 0{f 



{N + sy 9 

Defining = Z4/ZG and using (jS^) , we obtain 
loM = Pig) — — =9 



dg 



25515 (AT + 8)2 



1 + ^ X 0.045007 + 0(^2) 



(163) 



The critical exponent a is obtained by evaluating 704(^7) at the fixed-point value of the 
coupling, i.e. 



= ioAq*)- 



(164) 



The scaling function ^4(?/) is obtained from the zero-momentum renormalized function 
Toifi- Indeed setting 



04,-R 



Qi{k)fi{g,y), 



we have 



9A{y) = lim U{g,y). 



9^9 



(165) 



(166) 



The expansion of f4{g, y) is 
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3 iV + 2 SOtt^ f 



[Q4 J2ik, 1) - g4 J2ik, l)U=o] 



(167) 



+^' (Ar + 8p"9' ^^/^^^ " ^^-^s.il^, 1) - 167rJ3,2(^, 1)J - = 0)| + 0(^^), 

By expanding f^lg, y) in powers of y around y = 0, one finds 



, 2 N + 2 , 
di = g 



di 



and 



di 
d2 
ds 



380 



168399 
3076 

19702683 
3112 



253320210 



1 + ^ X 0.105400 + 0(^2) 
"l - ^ X 0.355629 + 0(^2) 
1 - ^ X 0.696450 + 0{f) 



etc. 



(168) 



(169) 



In order to get estimates of a and of the coefficients di from the corresponding series, 
we have employed the procedure outhned in Section [111 C| . Results for a are reported in 
Table ^IT|, and for di in Table ^. 



E. An e-expansion analysis 

The exponents ri2i introduced in Sec. [11 A| can be evaluated in the e-expansion of the 
corresponding 0^ theory. A rather simple two-loop calculation gives 

_ 1 A^ + 2 2 , r^r 3 
''~2(iV + 

and 



+ O(e^) (170) 



We recall that the exponent i] is known to O(e^) 0. For a one then finds 

7 N -\-2 

a = ri-V^ = - / ,, e2 + 0{e'). (172) 
20 (AT + 8)2 ^ ' ^ ^ 

To compute g^iy) , we consider the renormalized two-point one-particle irreducible func- 
tion with an insertion of the operator 04(x), see (^9]). To order 0{e^) we find 

The function Js(/i;, m) is the finite part of the integral 
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with the MS prescription. The expansion of ^4(2/) in powers of y gives 



2 ^^ + 2 



and 



^^ = 'jf;^2d^ + 0{e^) (175) 



= -0.00354500..., 
d2 = 0.00011715..., 

4 = -0.00000599..., (176) 



etc... 



F. A strong-coupling analysis 

Anisotropy in the two-point function can be studied at finite by analyzing the strong- 
couphng expansion of its lowest non-spherical moments. 

In order to compute a, the correction to the Gaussian value of p, we analyze the strong- 
coupling expansion of the ratio g4,o/''^2, which behaves as 

^ ~ ~ (T - Te)"^ (177) 

for T Tc. We recall that in the 1/A^ expansion u = 1 + 0(1/A^), and for = 0, 1, 2, 3 
u ~ 0.588, u ^ 0.630, u ~ 0.670, u ^ 0.705 respectively 0. DPA's and lA's of the available 
strong-coupling series of the ratio g4,o/'^2 on both cubic and diamond lattice turned out not 
to be sufficiently stable to provide satisfactory estimates of a at any finite value of N. 

A more satisfactory analysis has been obtained by employing the so-called critical point 
renormalization method (CPRM) |^8|. The idea of the CPRM is that from two series D[x) 
and E{x) which are singular at the same point Xq 

D{x) = J2 dix' ~ (xo - x)-\ 

i 

E{x) = J2e,x' ^{xo-x)-% (178) 



one can construct a new series by 



F{x) = J2-x\ (179) 



The function F{x) is singular at x = 1 and for x — > 1 behaves as 

F(x) ~ (1 -x)-"^. (180) 
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The exponent is related to S and e by 



0=1 + 5-6. (181) 

Therefore the analysis of F{x) provides an unbiased estimate of the difference between the 
critical exponents of the two functions D{x) and E{x). Moreover the series F{x) may be 
analyzed by employing biased approximants with a singularity at Xc = 1. 

By applying the CPRM to the strong-coupling series of ^4^0 and m2, one can extract an 
unbiased estimate of cr, by computing the exponent (p = 1 — au from the resulting series at 
the singularity xq = 1. We analyzed this series by biased lA's, considering those indicated in 
Eq. ( p,18| ). The estimates of a we obtained confirm universality between cubic and diamond 
lattice, although the analysis on the diamond lattice led in general to less stable results. 
In Table |V111], for selected values of N, we report our estimates of a, which are essentially 



obtained from the analysis on the cubic lattice. In order to derive a from au, which is the 
quantity derived from the strong-coupling analysis, we have used the values of u available in 



the literature. See e.g. for an updated collection of results obtained by various numerical 
and analytic methods. The errors we report are rough estimates of the uncertainty obtained 
by considering the spread of all the analyses we performed. The values of a are very small 
at all values of A^, and at large N, say ^ 10, they are consistent with the corresponding 



0(1 /N) prediction, cfr. Eq. (154). By analyzing separately the strong-coupling series of q^^ 
and m2 in the case of a, and taking the difference of their exponents, one obtains consistent 
but less precise results. 

In order to estimate the first non-trivial coefficient di of the expansion of g4,{y), see 
Eq. ( |121| ), one may consider the quantity ri defined in Eq. ( |1 23| ) . However, as we did for 



the analysis of q in Sec. |Tl|, it is better to consider another quantity ri which is defined so 
that ri = for N = 00 for all j3 < j3c- For the cubic lattice 

fi = 2 - ^ ■ ^ + 182 

while for the diamond lattice 

In the critical limit fi rfi. On the cubic lattice the available series of ri has the form 
J2t=o '^iP^ ^ except for = 1 where it has the form J2i=o '^iP'' ■ These series can be 
derived from the strong-coupling expansion of G{x) presented in App. |^. On the diamond 
lattice the available series of ri has the form I]i=o ^«/^*' except for A^ = 1 where it has 
the form P^ YllLo diP''- The results of the analysis are reported in Table Universality 
between the cubic and diamond lattice is again substantially verified, although the diamond 
lattice provides in most cases less precise results. The value of di is very small for all A^. At 
large- A^ the strong-coupling estimate of di is in good agreement with the large- A^ prediction 
( |155|) . The estimates are also in satisfactory agreement with the results obtained from the 
^f-expansion and the e-expansion. 

We have also obtained estimates of = rj — rjQ, i.e. the anomalous dimension of the 
irrelevant operator 
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Oti{x) = 0(ifl{x) = s{x) ■ Qe{d)s{x), 



(184) 



by analyzing the strong-coupling series of g6,o- For the cubic lattice the exponent erg is 
determined from the critical behavior of the ratio qe.o/'m2, since 

^ ~ ~ (T - T.y'", (185) 



Notice that Eq. (|185|) is not valid on the diamond lattice, since here gg.o receives contributions 



from two operators, from Oq{x) and from the leading irrelevant operator responsible for the 
parity breaking. Results for various values of are reported in Table |X|. They were 
obtained by applying the CPRM to the series of ge.o and m2, and by analyzing the resulting 
series by biased lA's. Like a, ctq is small for all values of A^. At large A^ our estimates 
compare well with 

^6 = + O f 4 V (186) 



337r2A^ \m 

Finally we compute pp for the diamond lattice. For a Gaussian theory pp = 3 and thus 
93,0 ~^ constant for Mq —>■ 0. In general, at finite values of A^, we write Pp = 3 + cTp. The 
exponent Up is determined from the critical behavior of ^3 g 

q,o-M^r (187) 



In order to estimate Up, we applied the CPRM to the series q^^ and x- We found < o"p 
0.01 for all A^ > 0. 



< 



G. The two-dimensional Ising model 

We conclude this section by considering the two-dimensional Ising model, for which 
we present an argument showing that the anomalous dimension of the irrelevant operators 
breaking rotational invariance is zero. 

Let us consider first the square lattice. In this case, for sufficiently large values of |x| (in 
units of the lattice spacing) the asymptotic behavior of G{x) on the square lattice can be 
written in the form 



where p = X]/i4sin (p^/2 



.2^2 ,.1/4(1 + ^2)^2 



and 



Z{(3)= (l-zy-Az' ' ^ ' ^ , (189) 



(1 + ^2)2 

M2(/3) = l-L_^-4, (190) 
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and we have introduced the auxihary variable 



z{p) = tanhp. (191) 

This shows that at large distances the breaking of rotational invariance is identical to that 
of the massive Gaussian model with nearest-neighbor interactions. Therefore p = 2 exactly. 

This value of p is confirmed by a strong-coupling analysis of the moments q^^rn using the 
available 21st order strong-coupling series [^]. In particular, on the square lattice we found 

lim ^ = 1 (192) 

within an uncertainty of O(10~^). 

A formula analogous to Eq. (|188|) has been conjectured in Ref. BH] for the Ising model 



on triangular and honeycomb lattices. Thus, also on these lattices, the pattern of breaking 
of rotation invariance (and parity in the case of the honeycomb lattice) should be that of the 
corresponding Gaussian theories, which have been described in Sec. [IV B| . If the conjecture 



of Ref. ||36| is correct, we have p = 4 for the triangular lattice and Pp = 3 for the honeycomb 
lattice. 

Again, an analysis of the strong-coupling expansion of G{x) on the triangular and hon- 
eycomb lattices supports convincingly this conjecture. 



H. Summary 

For lattice models with 0{N) symmetry we studied the problem of the recovery of 
rotational invariance in the critical limit. Anisotropic effects vanish as Mq, when Mq 0. 
The universal critical exponent p, which is related to the critical dimension of the leading 
operator breaking rotational invariance, turns out to be 2 with very small A^-dependent 
corrections for the lattices with cubic symmetry. Notice that this behavior is universal and 
thus should appear in all physical systems which have cubic symmetry. The reader should 
note that p is different from the exponent uj, which parametrizes the leading correction to 
scaling and which is related to a different rotationally-invariant irrelevant operator. Models 
defined on lattices with basis, such as the diamond lattice show also a breaking of the parity 
symmetry. We find that these effects vanish as Mq , with Pp ^ 3 for all values of A^. 

We have also calculated the universal function g4{y). For ?/ < 1, we find g4,{y) ~ 1 with 
very small corrections. 

In our study we considered several approaches, based on 1/A^, g-, e-, and strong-coupling 
expansions. All results are in good agreement. 

In two dimensions we showed that p = 2 for the square lattice for all A^ > 3 and A^ = 1. 
We conjecture that this is a general result, valid for all values of A^. Similar arguments apply 
to the triangular (honeycomb) lattice: we conjecture p = 4 (resp. pp = 3) for all A^. 
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APPENDIX A: THE LARGE- AT LIMIT ON THE CUBIC, F.C.C. AND 

DIAMOND LATTICE. 



In this appendix we present the large- hmit of the two-point function in O(A^) a model 
with nearest-neighbor interaction on cubic, f.c.c. and diamond lattices. 

1. The cubic lattice 



The large- two-point Green's function on the cubic lattice is ||6y] 



(Ai) 

where 

i 

^ = 2{l-cosk,), (A2) 
and A^s is the number of sites. In Eq. ( [Al[ ) z = Mq, where Mq is the zero-momentum mass 



scale. The lowest spherical and non-spherical moments of G{x) are reported in Table ^1. 

The relation between (5 and z is determined by the condition G{0) = 1. In the infinite 
volume limit one has 

r d^k 1 , . , 



An expression in terms of elliptic integrals can be found in |^]. The critical point /?c is 
obtained by setting 2; = 0. One finds 



2 2 

= (18 + 12V2 - 10V3 - rVe) a: (2 - 2 V2 - 3 + Ve) =0.252731..., (A4) 



7r2 



where K is the complete elliptic integral of the first kind. The expression for (3c in terms of 
F-functions reported in (and quoted in [Q) is unfortunately wrong. 

Concerning tha mass-gap estimator introduced in Eq. ( |115| ), it is easy to check that for 
A^ = oo and for 13 < 13^ 



4 = ^- (*^' 



independently of (3. 
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2. The face-centered cubic lattice 



The f.c.c. lattice is a regular lattice, thus the solution of the corresponding Gaussian 
model does not present difficulties. Its coordination number is c = 12 and the volume per 
site is = (in unity of a^, where a is the length of a link). 

The sites x of a finite periodic f.c.c. lattice can be represented in cartesian coordinates 

by 

X = hffi + + ^3^3, 

r/i = ^(0,1,1), 

7/2 = ^(1,0,1), 

7/3 = ^(1,1,0). (A6) 

It is easy to derive the massive Gaussian propagator and therefore the large- two-point 
Green's function, which is 

where 

A{k) = 2 (3 - cos ^ cos ^ - cos ^ cos ^ - cos ^ cos ^) , (A8) 

and Ng = L1L2L3 is the number of sites. The momenta k run over the reciprocal lattice: 

- 27r ^ 27r ^ 27r ^ 
k = —mipi + —m2P2 + -;-m3p3, 
Li L2 

TTli = 1, ...Li, 

Pl = ^ (-1,1,1), 

P2 = ^ (1,-1,1), 

P3 = ^ (1,1,-1). (A9) 



Again z = Mq. The lowest moments derived from Eq. ([A7|) are reported in Table [VI 



The relation between (3 and z is determined by the condition G'(O) = 1. In the infinite 
volume limit one obtains 

/?=ir^ \ , (Aio) 

^ 2i-^ (27r)3 2D(g) + z' ^ ^ 

where 

D{q) = sin^ f + sin^ f + sin^ f + sin^ + gin^ 3a + gin^ «. (All) 
By setting z = one finds ^2 



39 



/?.= ^^ = 0.112055.... (A12) 

It is worth noticing that the function eo(A;^), defined in Eq. (p^, turns out to be the 
same for cubic and f.c.c. lattices, as a consequence of some trivial symmetries of the angular 
integration. Therefore the first few spherical moments are equal as functions of Mq even 
off-criticality, as shown in Table [V^ . 



3. The diamond lattice 

The diamond lattice has coordination number c = 4. It is not a regular lattice, because 
not all the sites are related by a translation. It consists of two interpenetrating f.c.c. lattices, 
and can be regarded as a f.c.c. lattice with a two-point basis. The absence of translation 
invariance causes a few subtleties in the analysis of models defined on it. One cannot define a 
Fourier transform which diagonalizes the corresponding Gaussian propagator. Neverthless, 
observing that sites at even distance in the number of lattice links form regular f.c.c. lat- 
tices, one can define a Fourier-like transformation that partially diagonalizes the Gaussian 
propagator (up to 2 x 2 matrices). 

Setting the lattice spacing (i.e. the length of a link) a = 1, the sites if of a finite periodic 
diamond lattice can be represented in cartesian coordinates by 

X = x' + pfjp 

x' = hVi + kV2 + km, 

/i = l,...Li, p = 0,l, (A13) 

where 

r/p = ^(l,l,l), 
r/i = ^(0,1,1), 

7/2 = ^(1,0,1), 

7/3 = ^(1,1,0). (A14) 

The total number of sites on the diamond lattice is A^^ = 2L1L2L3, and the volume per site is 
Vs = (in unit of a^). The variable p can be interpreted as the parity of the corresponding 
lattice site: sites with the same parity are connected by an even number of links. 

The two sublattices identified by x+{li,l2,h) = x{li,l2,h,0) and X-{li,l2,h) = 
x{h,l2,h,^) form two f.c.c. lattices having A^^ = A^s/2 sites. Each link of the diamond 
lattice connects sites belonging to different sublattices. It is convenient to rewrite a field 
0(x) = (j){li,l2,h,p) in terms of two new fields 0+(x+) = 0(a^+) and 0„(x_) = 0(x„) de- 
fined respectively on the sublattices x+ and x_. A finite lattice Fourier transform can be 
consistently defined by 

x± 

MS±) = ^Y.^^''''^Mk), (A15) 

" k 
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and the set of momenta is 



27r ^ 

rrii = 1, ...Li, (A16) 



where 



ZTT 




ZTT ^ 




TJliPi 


-f —m2p2 

I JO 


= 1 


...Li, 




pi 


= vi 

4 


(-1,1,1), 


P2 


= 

4 


(1,-1,1), 


P3 


= f 


(1,1,-1), 



(A17) 
so that 

- ^, 27r , 27r , 27r , , 
K ■ x = —mill + —7712/2 + -j-m^h- (A18) 
Li L2 -L3 

The large- iV hmit of the two-point function is the massive Gaussian propagator defined 
on the same lattice. By using the Fourier transform ( [A16|) , straightforward calculations allow 
to derive a rather simple expression of the massive Gaussian propagator, and therefore of 
the large- two-point function: 

G{x;y) = G{x',p^;y',py) = 



where 



(A20) 



^ k- ^ k- 
H{k) = n cos ^ - 2 n sin 

i=i VO i=i VO 



and 

/\{k) = ?>{l-\H{k)\ 

= i (3 - ^ % - % * - * • (A21) 

Notice that A(A;) has the same structure of the inverse propagator of the f.c.c. lattice, cfr. 
Eq. (|A8|). One can easily verify that z = Mq, where Mq is the second-moment mass. Using 
Eq. ([A19|) one can derive the expression of the lowest moments of G{x) reported in Table [VI . 

The relation between (3 and z is determined by the condition G{0) = 1. In the infinite 
volume limit one can write 



4 7-. i27ry ^^D{q) + z {1 + 1 



12^ 



where D{q) has been already defined in Eq. (|A11|) . Comparing Eqs. ([A22|) and ( |A1CI|) for 
z = 0, one notes that /3f ^'""^'^ = Ap^-^-^-. Therefore for the diamond lattice (3c = 0.448220.... 
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In order to define a mass-gap estimator on the diamond lattice, one may consider the 
wall- wall correlation function defined constructing walls orthogonal to w = -^(—1,1,0), 
which is the direction orthogonal to two among the links starting from a site. We define 

G^{t = x-w)= J2 <^(^) (A23) 

t=cst 

where the sum is performed over all sites with the same t = x ■ w = -^{h — h) (the 
coordinates of the sites x are given in Eq. ( [A13| )). Using Eq. ([A19| ) one can easily prove 
that G'^(t) enjoys the property of exponentiation. The mass-gap /i can be extracted from 
the long-distance behavior of Gw{t). For t ^ 1 

Gy,{t) oc e"^* (A24) 

In view of a strong-coupling analysis, it is convenient to use the following quantity 

4 / I- 



Mi = - (^cosh^f /i - 1 j (A25) 

which has the property /i for /i — and has a regular strong-coupling series. In the 

large- iV limit and for (3 < (3c 

Ml , , , 



APPENDIX B: PERTURBATIVE EXPANSION OF SCALING FUNCTIONS 



In this appendix we present a simple derivation of all the results that are needed in 
order to construct explicitly the 1/iV, and g- and e-expansions up to three loops presented 
in Sec. 



[II 



Our starting point is the observation that most of the two- and three-loop 
calculations needed in the relevant perturbative calculations are included, apart from rather 
trivial algebraic dependences on A^, in the one- loop calculation of the 1/A^ expansion for the 
two-point function. As we shall show, the 1/A^ results can be expanded in g and e in order to 
recover all the corresponding contributions. Let's therefore start with the evaluation of the 
renormalized self-energy to 0{1/N) in arbitrary dimension d and for arbitrary bare coupling 
go in the A^-component 0^ theory. 

We introduce the dressed composite propagator (geometric sum of bubble insertions in 
the 0^ vertex): 



1 r d^p I 1 3 

2 J {2TiYp^ + m? {p + kf + m? ^ Wg 



m'-'^ ^ /\;\y) + (Bl) 



where y = k'^/m?, and we have defined the (zero- momentum subtracted) dimensionless 
renormalized dressed (inverse) propagator: 



-\y) = A'\y,go)-/\-\Q,go 
1 r d'^p 1 
~2J {2'KYp'^ + m? 



{p + kY + m? + m? 



(B2) 
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and the four-point (large- A^) coupling renormalized at zero momentum 



(B3) 



where we have rescaled the coupling for convenience, generalizing a rather standard three- 



dimensional prescription. The integration (|B2|) can be explicitly performed, and one obtains 



2 (4ir)*2 



F 2 



2' 2' 2' y + 4 



2(47r)^/2 



(B4) 



which is a regular function of d for all < 4. 

The renormalized contribution to the self-energy can now be evaluated by the formal 
expression 



d 



(B5) 



(T{y,g) = m 



2_rf 2(47r)''/2 ^ 



r 2 



{2ttY 1 + g5r{p'^/m?) {p + kf + m? 



(B6) 



and the subtractions that are symbolically indicated in Eq. ( P5| ) must be done before per- 
forming the integration in Eq. ( P6|) in order to obtain finite quantities in all steps of the 
derivation. To this purpose, it is convenient to perform first the angular integration, by 
noticing that 



2-d 



m- 



;f{p'/m') = 2B (f , 2 - f ) {zf^~' dzf{z)h{z, y) 

(B7) 



where 



h{z,y) 



d9- 



sm( 



-.d-l 



b[^,\)Jo z + y + l + 2^cose' 
The subtraction indicated in Eq. (|B5| ) now simply amounts to replacing in Eq. ( p6| 



[B8) 



d 

h{^, y) — ' h{^, y) - h{z, 0) - y) \y=o 



h{z,y) 



+ 



y 



Azy 



i + z ' {1 + zy d{i + zf' 

The relevant results for some integer values of d are: 

4 



(B9) 



Sriy) 

h{z,y) 



- 1, 



(4 + y) 

z + y + 1 

{z + y + iy - Azy 



(BIO) 
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for d = 1; 



h{z,y) = \iz + y + lY-Azy 



-1/2 



(Bll) 



for d = 2] here ^ 



1 + 



h{z,y) 



2 ^yy 

— — arctan-^^ 1, 

1 ^^z + y + 1 + 2^ 



for d = 3: 



limr(2-f)5,.(y) 



lim/i(2;,?/) 

d— >4 



z + y+1-2^ 



2; + !/ + 1 - 



(B12) 



2zy 



z + y 



Azy 



(B13) 



Eqs. (p5|) and ( pq) are now ready for our purposes. Let us immediately notice that, 
since we know the leading large- iV fixed point value of g, which happens to coincide with 
the Ng^ —>■ oo limit of Eq. ([B3|), we may replace 
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9 



(B14) 



and find the contribution to the scaling function go{y), which in turn is simply the 
continuum non-linear a model evaluation of the self-energy. This is the way Eq. (|68D is 
generated, by setting o? = 3 in the general expression. Eq. ( [B13| ) shows that, at least in the 
naive — > 4 limit, there is no 0(l/A^)-non-Gaussian contribution to the self-energy scaling 
function. 

Eq. (|E6| ) is also the starting point for the g- and e-expansion up to three loops. It is 
indeed straightforward to obtain a representation of the leading 0{1/N) contributions to 
the self-energy as a power series in g: 



My, 9) = -fMy) + fvM + 0{f) 

where we have defined the functions 



(B15) 



5 f,2-fWo 



oo d 

Z'2~ dz[5r{z)] 



n-1 



h{z,y) 



1 



+ 



Azy 



l + z il + z) 



d{l + zf_ 

(B16) 



and we exploited the trivial consequence of the definition Eq. ( pi6|) : ^i{y) = 0. Restoring 
the correct dependence on for arbitrary (and not only very large) values of in front of 
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the functions and ^93 is now simply a combinatorial problem, whose solution leads to the 
complete three-loops result 



N + 2 
(iV + 8) 



^ , , _3 (iV + 2) 



(B17) 



We must keep in mind that the functions ^Pn{y) carry a dependence on the dimensionality 
and the scaling function go{y) is the value taken by /(y, g) when evaluated at the fixed point 
value g* of the renormalized coupling, where g* is in turn a function of the dimensionality 
and it is obtained by evaluating the zero of the /5-function. We may now choose two different 
strategies. The first simply amounts to fixing d to the physical value we are interested in 
and replacing g* with the numerical value (possibly evaluated by a higher-order expansion of 
the /9-function at fixed dimension). We may however decide to expand the functions ^Pn{y) 
in the parameter e = 4 — d around their value at d = 4, perform a similar expansion for the 
/9-function, then finding g* as a series in e 



r(d/2)(47r)'^/' 




3(3A^ + 14)^ 
(AT + 8)2 



+ 0(6^). 



(BIS) 



The functions ip2{y) and ip3{y) we have introduced in Sec. [Ill C| are strictly related to ip2{y) 
and (p3{y) calculated for d = 3, indeed 



¥^2{y) = V2{y)\d=3 

v^{y) = [^3(1/) - 2^2{y)]d=3 ■ 



(B19) 



APPENDIX C: 0{1/N) CALCULATION OF g^iy) 

In this appendix we present the calculation to order of the scaling function gi{y). 
We will first study the small-p behaviour of the auxiliary propagator 

1 , I r d'^q 1 

We are interested in the behaviour for p 0, Mq with arbitrary ratio p/Mc- For 
d < 4, the only case we will consider, the leading order is simply given by the continuum 
expression 

A -If' ^ f '^'^^ ^ 
° ^^^^2 leant {2'kY [{q + pf + Ml] (g2 + ) " ^^^^ 

Here and in the following, when we will append the subscript "cont" to the integral, we 
will mean that the integration domain is the whole (i-dimensional space; otherwise the 
integration should be intended over the first Brillouin zone. We want now to compute the 
first correction. Using 

f = q^ + eoM^f + e4,oQ4(g) + O(g^), (C3) 
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we could expand ( |C1D in powers of q. This however cannot be the correct answer as the 
integral one obtains in this way is ultraviolet-divergent. There is however a standard way 
out. Introduce a fixed cut-off A and define a sharp-momentum regularized quantity 



J(g,p;MG,A) 



((g + p)2 + M2)(g2 + M2)2 
+ Ml- 2{eo,2iqy + e,,oQ,{q)) + 2eo,2(g' + M'^feHq + pf - A'] 



(C4) 



where 9{x) is the step-function [0{x) = 1 for x > 0, 6'(x) = for x < 0]. Then we rewrite 
JTD as 



Z, J a 



cfiq 

cont (27rp 
1 



XlatM) 



+ - 



d'^q 

cont 



I{q,p;MG,A) 



jI{q,p-MG,A) 



(C5) 



where Xiattil) is the characteristic function of the first Brillouin zone, i.e. Xiattil) = 1 if 
q G [— vr, ttY, Xiattio.) = otherwise. For p — and Mq (but only for 2 < c? < 4), we 
can simply set Mq = p = in the first integral thus obtaining a constant which depends on 
A and on the specific lattice Hamiltonian. We must now compute the second integral. A 
completely straightforward calculation gives 



cont 



d'^q 

W) 



Jiq,p;MG,A) = A,\p) -eo,2 



d'^q 



-64,0- 



2<A2 {27rY g2 + ^2 
A{p') + Bip')A,\p) 



- AMlA~\p) - Uf-^ 



dp'^ 



(C6) 



where 



(47r)-'^/2r(3 - d/2) p\d + 8) + 36MI 



B{p' 



2(d- l)(d-4) 
1 1 



p2 + 4M2 



G 5 



(4 - d){d+2){py - 8{d - lOyM^ + 144M^ 



8(d- I)p2 + 4M2 
Summing the two terms we can rewrite the auxiliary propagator as 

A~\p) ^ A^\p) + C{Mg) + eo,2/o(p', Mg) + e^^^Q ,{p) Up\ Mg) 



where 



C{Mg) 
Up\Mg) 

Mp^Mg) 



d'^q 

cont {2ttY 



-M^4A^\p) + M'a 



Xiattjq) ^ 

dA^\py 



2eo,2 I -T 



g2 q2 + Ml)_ 



Qp2 



{p 



\-^{A{p') + B{p')A-,\p)). 



(C7) 
(C8) 

(C9) 

(CIO) 
(Cll) 
(C12) 
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The functions foip^^Mc) and fiip^^Mc) are universal: they do not depend on the lattice 
structure. The lattice dependence is contained in the two constants eo,2 and 64^0 and in 
C{Mg) which depends on the explicit form of the lattice hamiltonian. 

Let us now consider G~^{k, Mq) to order We will first compute Z^{Mg) = 

(74(0, Mg) which can be easily obtained from 



Z4(M( 



d + 2 



G 



-Qid/dk) G-'ik)\k=o 



2Ad{d - 1) 

Using the expression for G~^{k) to order cfr. formula (|137|) , we get 



Z^{Mg) = 64,0 + 



1 d + 2 r d'^q 
N24:d{d-1) J (2^ 



(C13) 



(CM) 



where 



6z4{q) 



D{q? 
8 

36 
6 



d 1 ^--v 3 ^--v 

1+2 ^ '^''^^^ ~ d + 2 ^ 



d — 1 ^ 3 

2^ W,,,,,,.Wn. - 



d + 2^ ^ d + 2 



D{qf 
72 



d-l 
d+2 



d-l 



1 



E 



w,. 



d + 2^ d + 2 , 



Y,{Wf,^w^^ + 2wl^) 



24 



D{qy 



d- 1 



E 



d+2 



E 



2 2 



(C15) 



Here D{q) =f + M^, = d^^f 



dudaq^ and so on. From this expression one can 



easily compute the exponent 774: indeed one must expand 6z4,{q) and A(g) for g — > and 
keep only those terms that behave as Ao(g)((3'^ + M^)~^. We obtain in this way the results 
iW WM for / = 2. 

We want now to compute the scaling function g^iy)- First of all notice that 



gA{y,MG) 



Ma 



d'n{k)^G~\k,MG) 



{k 



(CI6) 



where d'^Q{k) indicates the normalized measure on the {d — l)-dimensional sphere and we 
obtained 



d'^Q{k)QA{kf 
J d''n{k)Qeik) 



24(d- 1) 



(rf + 2)2(d + 4)(d + 6) 
2_ 720{d - 2) {d - 1 



{kY=MikY 



{d + 2){d + 4)2(rf + Q){d + 8)2(d + 10) 



(kY^^f6ik 



2\6 



(C17) 
(CIS) 
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(We shall not need Eq. ( |C18|) , but we quote it for further reference.) Using then Eq. (|137|) 
we get 



NNaJ {k^r J {^t^Y {q + k) + Ml 



(C19) 



so that 



My) = l + ^^-. hm let'em 



2\4 



(2vr) 



A(g) 



d + 2 



{q + k) + Ml 24rf(rf - 1 



-Q^{k)^{q)5z^{q) 



(C20) 



Now because of the subtraction, we can take the limit Mq — > by expanding the integrand 
in powers of q and using ( |C9|) . We now integrate over the angular variables. In particular 
we need to compute the following two integrals: 



Using the techniques of Ref. [T^l we get 

h{k\q\Ml) = {k'q'fl'F,,,{z) 



(C21) 
(C22) 

(C23) 



h{k\q\Ml) = 
where we have defined: 



2g2 



/ 2 \ 1/2 ^2 



(C24) 



g2 + fc2 + 



G 



2^q^ 
r((rf-l)/2)(d + /-3)! 



(C25) 



.l)^e-('^-^W^(.^-l)('^-^)/^g!t^^/3^)/,(z) . (C26) 



Here Q'^{z) is the associated Legendre function of the second kind (see Ref. ||6^, Sec. 8.7 
and 8.8). Notice that for / = the functions F^fi^z) are related to h{z,y) defined in ( |B8D . 
Indeed 



q^k"^ h 



Ml' Ml^ 



(C27) 



In two dimensions it is easy to see that 



48 



For d — 3 we have 



(C28) 



-Qi{z), 



(C29) 



and ior d — A 



— 1 — z 



(C30) 



Putting all together, we get finally 



11 

9i{y) = 1 - T7-I 



Ny^J (27r) 



Ao(g)V4(g', 1) + Ao(g)/2(y, g', 1) - subtr (C3i) 



where "subtr" indicates the integrand computed for y ^ 0. As expected the final result is 
universal: any dependence on the lattice hamiltonian has disappeared. 



APPENDIX D: STRONG-COUPLING EXPANSION OF G{x) ON THE CUBIC 

LATTICE 

Presenting l-th order strong- coupling results for the two-point Green's function would 
naively imply writing down as many coefficients as the number of lattice sites that can be 
reached by a /-step random walk starting from the origin (up to discrete lattice symmetries). 
It is interesting to notice the relationship existing between the number rii of lattice points 
(not related by a lattice symmetry) that lie at a given lattice distance / from the origin 
and the number of independent lattice-symmetric functions Q2m{k){k^y~"^- One can easily 
get convinced that, on a hypercubic lattice, the number of functions Q2m(^)(^^)' ™ is the 
same as the number of monomials of total degree / in the d variables that are not related 
by a lattice symmetry (that is, the number of independent, homogeneous lattice-symmetric 
degree-Z polynomials in the kf). This number in turn is equal to that of the partitions of 
I into d ordered non-negative integers, and this is nothing but the number of independent 
lattice points at a lattice distance I (where ordering insures independence by elimination of 
copies obtained by permutation). As a corollary, the relationship pi = rii — ni-i holds for 
arbitrary d on hypercubic lattices. A generating function for ni, for a given value of d, is 

oo d 1 

1=0 n=l ^ 

implying the asymptotic behavior 

In the case of three-dimensional hypercubic lattices, one can show that pi = |_//6J + 1 
with the exception of Z = 6A; -|- 1 in which case pi — k, while ni is the integer nearest to 
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(/ + 3)^/12 and the sum X^Leven'^i the integer nearest to (/ + 4)^/72. This would mean 
roughly (/ + 4)^/72 coefficients for the l-th order of the strong-coupling expansion on the 
cubic lattice. This number can be sensibly reduced (asymptotically by a factor 27 on the 
cubic lattice), without losing any physical information, by noticing that the inverse two-point 
function, when represented in coordinate space, involves only points that can be reached by 
a [//3j-step random walk. This fact can be traced to the one-particle irreducible nature of 
the inverse correlation. As a matter of fact, instead of the 93 coefficients needed to represent 
the 15-th order contributions to G{x), only 8 coefficients are enough for the corresponding 
contribution to the inverse function G~^{x), which we construct by the following procedure 
(a similar representation was used for the Ising model in a magnetic field in |^). 

We introduce the equivalence classes of lattice sites under symmetry transformations and 
choose a representative y for each class: whenever x ^ y then G{x) = G{y). We define the 
"form factor" of the equivalence class 

Q(?/) = Ee'"^ (D3) 

and represent the Fourier transform of G{x) according to 

Gip) = j:Qiy)G{y). (D4) 

y 

The momentum-space inverse Green's function is defined by 

G-\p)G{p) = 1. (D5) 

Therefore its inverse Fourier transform enjoys the symmetries of G{x) and satisfies the 
relationships 

G~\p) = Y.^'^^G-\x) =Y.Q{y)G~\y). (D6) 

X y 

In practice we exploit the property 

Qiv)Qiy) = J2niz;v,y)Qiz) (D7) 

z 

where 

n{z;v,y)= ^ S^^u+x (D8) 

are integer numbers, and reduce the problem of evaluating G^^{y) to that of solving the 
linear system of equations 

J2G-\v)M{v,z) = 5,,o (D9) 

V 

where 

M{v,z) = J2G{y)n{z;v,y). (DIO) 
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When expanding in powers of (3, the system takes a triangular structure and , as expected, 
it admits a solution whose non-trivial terms are only those corresponding to the equivalence 
classes of sites that can be reached by //3 random steps. 

Solutions for G~^{x) can be found for arbitrary A^. In Table ^we only exhibit G~^{x) 
for = 0,1,2,3,4 and 16. We choose a representative of the equivalence class by the 
prescription Xi > X2 > X3 > 0. We may notice as a general feature that in the class 
represented by xi > 1, X2 = X3 = the first non-trivial contribution is of order 3xi + 2 
{3xi + 4 when N = 1). When = 0, 1 a number of seemingly non-trivial coefficients turn 
out to be zero. 



APPENDIX E: STRONG-COUPLING SERIES OF x AND ma ON THE 

DIAMOND LATTICE 

On the diamond lattice we have calculated the strong-coupling expansion of G{x) up to 
21st order. In the character-like approach, the possibility of reaching larger orders than on 
the cubic lattice is related to the smaller coordination number. However longer series do 
not necessarily mean that more precise results can be obtained from their analysis. This 
is essentially related to the approach to the asymptotic regime of the correspoding strong 
coupling expansion, which is expected to occur later on lattices with smaller coordination 
number. The 21th order series on the diamond lattice provide estimates of the exponents 
7 and u which are, as we shall see for N = 1,2,3, substantially consistent with the results 
obtained by analyzing series on cubic-like lattices (see for example Ref. where series to 
0(/?^^) for the cubic and b.c.c. lattice have been presented and analyzed), but less precise. 

In this appendix we report the 21st order strong-coupling series of x and m2 calculated 
on the diamond lattice, for = 1, 2, 3. 27th order strong-coupling series for A^ = 0, i.e. for 



the self-avoiding walk model, can be found in Ref. ||42|| . 



1. N = 1 



^2 104/33 . 4328 /^^ 12128 711328/3^ 132452 /3« 

y = 1 + 4(3 + 12(3^ + — + 100/5^ + — + — + — + — + 

A T-/^-r/^^g^ ^15 15 315 21 

49894088 230044448 /^^^ 20986492048/3^1 11593048528/3^2 

2835 ^ 4725 ^ 155925 ^ 31185 ^ 
6239638466896/3^3 40044715794736 /3i^ 381115667726672 /3i5 907261838473556 /3I6 

6081075 ^ 14189175 ^ 49116375 ^ 42567525 ^ 
635228192216156408 (3^'^ 5223277546855685888 (3'^^ 815815904018756584288 (3^^ 

10854718875 ^ 32564156625 ^ 1856156927625 ^ 

744572898253973823856/3^0 642020997051581736673936/921 

— H — + 0(B^^). (El) 

618718975875 194896477400625 ^ ^ 

^2 488 /33 2048/3^ 38888/3^ 417664/3^ 10027936/3^ 33306368 /3S 
m2 = 4 /? + 32 /?2 + '— + ^ + ^ + ^ + '— + ^ + 

^ 3 3 15 45 315 315 
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971601608/^9 15453950464 532482065296 /^^^ 4939730085376/5^2 

2835 ^ 14175 ^ 155925 ^ 467775 ^ 
196443743845456 /?i3 4168605624019328 /^^^ 188065240470724112 /?^5 

6081075 ^ 42567525 ^ 638512875 ^ 

561744708980235008 28352355075085440248 /3^'^ 

638512875 ^ 10854718875 ^ 

57966531061027107328 /3i^ 42081281751167641189216 /J^^ 

7514805375 ^ 1856156927625 ^ 

615759333006052918694656/^20 37696556941296724618984336/^21 ^.^22n 

9280784638125 ^ 194896477400625 ^ 

We have analyzed the series of x by using the [m/l/k] first order lA's with 

m + / + A; + 2 = 21, 

[{n - 2)/3j - 2 < m, /, A; < [(n - 2)/3] + 2. (E3) 

We have obtained Pc = 0.3697(1) and 7 = 1.238(14). An estimate of 7 can be also obtained 
by applying CPRM to the series cind as explained in Sec. |IVF| . By employing biased 
lA's (those indicated in ([E^), one finds 7 = 1.253(4). By applying CPRM to the series m2 
and and using biased lA's, one finds u = 0.645(4). These values of 7 and u are slithly 
larger than the available estimates obtained by other tecniques (field-theoretical approaches 
give 7 ~ 1.240 and z/ ~ 0.630), or strong coupling expansion on other lattices, but we would 
not consider them inconsistent. One should not forget that the error displayed does not 
take into account possible systematic errors (due, for example, to confluent singularities), 
but just the spread of the results of the various lA's indicated in (^3l). 



2. N 



889703 /^^ 2387483 /3i° 22968773 /J^^ 25617551/5^2 

— H — H — H — + 

60 60 216 90 

11516036093 /?i3 40849680041 /^^^ 520550507027 /J^^ 

15120 ^ 20160 ^ 96768 ^ 

3457894675397 /?i^ 495995794312009 /J^^ 

241920 ^ 13063680 ^ 

2188572410969059 (3^^ 173608313274399461 f3^^ 

21772800 ^ 653184000 ^ 

76543471229019871/320 5344313242348050991/321 

— + — + 0(3^^). E4 

108864000 2874009600 ^ ^ ^ 



„ ^4 7534/35 26488/3*^ 356305/3^ 289444 /38 
m2 = 4/3 + 32/32 + 162/3^ + 672/3^ + + + 777-^ + 5-^ + 
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18326503 /^^ 42659326 3125910649/5^1 1176454982/3^2 

— H — H \ h 

60 45 1080 135 

78423473449 /?i3 577822206313 /^^^ 108069034519903 /^^^ 

3024 ^ 7560 ^ 483840 ^ 

58770348791597 /3^^ 24384512261505001 /^^^ 

90720 ^ 13063680 ^ 

43660988509648999 /3^^ 9954936929950180901 f3^^ 

+ rTTTTTTTTTTrT^ — + 

q21 

22\ 



8164800 653184000 
1764942584095467281 /3^^ 1754883256361403082267 (3^^ 



40824000 ^ 14370048000 + '^('^ )• 

By performing an lA analysis of the series of x, one finds jSc = 0.3845(2) and 7 = 1.33(2). 
By applying CPRM to the series and and employing biased lA's, one finds 7 = 1.34(1). 
By applying CPRM to the series m2 and and using biased lA's, one finds u = 0.689(8). 
These results are substantially consistent with the available estimates of 7 obtained on other 



lattices and by other approaches (see e.g. Refs. 401 and 



3. N = 3 



^2 168 /53 468/5^ 9144 123456/5'' 65568/5^ 873708/5^ 
X = l + 48 + 12 8^ + — + — + — + — + — + — + 

128270568 /?9 11818853472 /J^o 2007117038928/5^1 5262987995856 /J^^ 

9625 ^ 336875 ^ 21896875 ^ 21896875 ^ 
1973906542032 25696714370736 277395071474138256 

3128125 ^ 15640625 ^ 65143203125 ^ 
5048136975344060076 /^^^ 1747312876419771883176 /^^^ 

456002421875 ^ 60648322109375 ^ 

35523883350405253078656 /^^^ 3207088211587054727672352 p^^ 

476522530859375 ^ 16678288580078125 ^ 

8294186293466843988864336/^20 285612671193686662161552/521 ^,^22n 

16678288580078125 ^ 221862716796875 ^ ^ 



^9 808 /53 3328/5^ 17240/5^ 1498496/5*5 4978592/5^ 15959296 /5S 

m2 = 4 /5 + 32 /52 + — + — + — + — + — + — + 

^^ 5 5 7 175 175 175 

391158744 /59 292871549952/5^0 8170771755824/5" 169326765636096/5^2 



+ + ^TTTTTTTT^ + ■ + 

+ 



1375 336875 3128125 21896875 

495146153921968 /?i3 7166586778308992 /5i^ 1747788082514945008 /5i^ 



21896875 109484375 9306171875 

243755148694999429888 p^*^ 91627122308762345759912 p^^ 

456002421875 ^ 60648322109375 ^ 

2022510678813614989101568 p^^ 197787508407138584345236512 p^^ 

476522530859375 ^ 16678288580078125 ^ 

21994072978677629556242688/520 34998691725014346631615751056 /521 Q(n22^ (-^^^ 
667131543203125 ^ 383600637341796875 ^ )• ( ) 
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By performing an lA analysis of the series of x, one finds j3c = 0.3951(2) and 7 = 1.42(2). 
We mention that singularities approximately as far to the origin as (3c have been detected by 
our analysis, indeed we found two singularities at /3 ~ ±i0.39. By applying CPRM to the 
series and x, and employing biased lA's, we obtained 7 = 1.42(1). By applying CPRM 
to the series 1712 and aiid using biased lA's, one finds u = 0.726(4). These results are 



slightly larger (and less precise) than the values obtained on other lattices (see e.g. Ref. |^ 
or by other techniques (see e.g. Ref. 0]), but substantially consistent. 
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TABLES 



TABLE L Numerical values of the coefficients /ij"'"'^ defined in Eq. (90) 



2 
3 

4 

5 



,(2) 



2 

'405 
4 



15309 
1 

59049 
4 

3247695 



(3) 



(4,1) 



(4,2) 



0.00949f25 
-0.000612784 

0.0000450060 
-0.0000035762 



0.000765804 
-0.00000341189 
-0.00000233206 

0.00000035769 



-0.0134856 
0.00102554 

-0.0000841861 
0.0000072651 



h 



(4,3) 



-0.0345992 
0.00253490 

-0.00020327 
0.00002390 



TABLE II. For several values of N and for the cubic and diamond lattice, we report the 
values of /3c we used in our strong-coupling calculations. We also report the fixed-point value g* 
of the rescaled zero-momentum four-point renormalized coupling, as obtained by field-theoretical 
methods. 



N 

IT 
1 

2 

3 

4 

8 
16 
32 
oo 



cubic 



diamond 



/3c 
/9c 
/3c 
/3c 
/3c 
/3c 
/3c 
/3c 
/3c 



=0.213492(1) 
=0.2216544(3). 
=0.22710(1) 
=0.231012(121 
=0.23398(2) [icr 
=0.24084(3) p 
=0.24587(6) 
=0.2491(1) 
=0.252731... 



/3c 

fie- 
/3c : 
/3c : 
/3c : 
/3c : 
/3c : 
P.- 
/3c : 



=0.34737(1) 

=0.3697(1) 

=0.3845(2) 

=0.3951(2) 

=0.4027(2) 

=0.4200(2) 

=0.4327(2) 

=0.4401(1) 

=0.448220... 



1.421(8) 
1.416(5) 
1.406(4) 
l-391(4iipi 
1.369 1221 
1.303 12^1 
1.207 iSI 
1.122 \22\ 
1 
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TABLE III. Estimates of C2, C3, and Sm — 1 for selected values of from various analyses of 
the strong-coupling scries on the cubic lattice (sec the text). An asterisk indicates that most of the 
approximants are defective, or, in the cases where numbers are not shown, that all approximants 



are defective, so that no estimate can be extracted. 





N 


PS 


PA 


DPA 


lA 


10* C2 





-0.6(1.8) 




* 


-1.0(3) 




1 


-3.0(4") 


-3.2f5) 


-3.0(6) 


-2.9(1) 




2 


-3.9(3) 


-4.1(7) 


-4.0(5) 


-3.9(1) 




3 


-4.1(1) 


-4.5(1.2) 


-4(3) 


* -3.7(4) 




4 


-4.06 


-4.4(6) 


-4.3(2) 


-4.09(3) 




8 


-3.4(1) 


-3.6(6) 


-3.59(9) 


-3.5(1) 




16 


-2.26 


-2.5(3) 


-2.4(1) 


-2.43(7) 




32 


-1.28 


-1.5(2) 


-1.46(9) 


-1.45(4) 


10*C3 





0.257 


0.12(2) 




0.125(3) 




1 


0.170 


11(41 


0.107(3) 


0.10(2) 




2 


0.205 


0.11(2) 


0.11(1) 


0.11(2) 




3 


0.212 


0.11(2) 


0.12(1) 


* 




4 


0.133 


0.12(1) 


0.12(1) 


* 




8 


-0.286 


0.07(4) 


0.10(1) 


* 




16 


-0.582 


0.04(3) 


0.070(5) 


* 




32 


-0.559 


0.02(2) 


0.041(3) 


*0.03(2) 


10*(5m-1) 





-2.03 


0(2) 


* 


* 




1 


-1.76 


-2(1) 


-3(1) 


* 




2 


-2.72 


-4(1) 


*-2(4) 


*-3.3(6) 




3 


-3.34 


-4(1) 


-4(3) 


-3.9(3) 




4 


-3.65 


-4.4(9) 


-4(3) 


-3.9(4) 




8 


-3.63 


-3.9(5) 


-4(2) 


-3.9(5) 




16 


-2.72 


-2.8(3) 


-3(1) 


-2.8(2) 




32 


-1.71 


*-1.8(3) 


-2.0(7) 


* 
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TABLE IV. Estimates of C2, C3 and Sm ~ 1 from various analyses of the strong-coupling 
series on the diamond lattice. An asterisk indicates tliat most of the approxiniants considered are 
defective, or, in the cases where numbers are not shown, that all approximants arc defective, so 
that no estimate can be extracted. 




N 


PA 


DPA 


lA 


10^02 





-2(1) 


* 


-0.6(3) 




1 


-3.0(3) 


-2.9(5) 


-3.0(2) 




2 


-4.1(4) 


-4.1(5) 


-4.36(4) 




3 


-4.6(3) 


-4.4(3) 


*-4.7(l) 




4 


-4.7(2) 


* -4.7(3) 


-4.8(2) 




8 


-4.0(2) 


* -3.9(2) 


-3.97(3) 




16 


-2.66(7) 


*-2.6(l) 


-2.65(2) 




32 


-1.52(3) 


*-1.51(7) 


-1.49(5) 


10*C3 





0.10(3) 


* 


0.099(6) 




1 


0.12(2) 


0.08(2) 


0.11(2) 




2 


0.12(4) 


0.08(1) 


0.12(2) 




3 


0.12(8) 


0.09(1) 


0.13(2) 




4 


0.1(2) 


0.09(1) 


0.12(2) 




8 


0.2(4) 


*0.08(2) 


0.1(1) 




16 


0.1(2) 


0.06(2) 


0.08(7) 




32 


0.0(1) 


0.04(1) 


* 


io-'CSm-i) 





0(1) 




* 




1 


-2.3(4) 


-2.2(3) 


-2.3(4) 




2 


-3.6(4) 


-3.4(2) 


-3.5(2) 




3 


-4.0(5) 


-3.9(3) 


*-4(l) 




4 


-4.3(7) 


-4.1(2) 


*-5(3) 




8 


-4.1(7) 


-3.6(3) 


-4.0(3) 




16 


-3.0(4) 


-2.4(2) 


-2.8(2) 




32 


-1.9(3) 


-1.5(2) 


-1.7(2) 
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TABLE V. Estimates of the coefficients C2 and C3, and of the mass-ratio Sm, from the analysis 
of the strong-couphng series of U2 and u^, M^/Mq (on the cubic lattice) and M^/AIq (on the 
diamond lattice). We report also results from the l/A'^-expansion, from the (^-expansion and the 
e-expansion. In the latter cases we give two numbers corresponding to the two choices: resumming 
R{x) or R{x)/x^. For N = 1 we also give the estimates from an "improved" resummation of the 



e-expansion which takes into account the exactly known results in two dimensions. 



N 






lO^Cg 


104(Sm - 1) 





cubic 


-1(1) 


0.12(1) 


0(2) 




diamond 


-1(1) 


0.10(1) 


0(1) 




g-expansion 


-3.29, -3.63 


0.108, 0.102 


-2.95, -3.50 




e-expansion 


-2.48, -4.26 


0.065, 0.114 


-2.55, -4.38 


1 


cubic 


-3.0(2) 


0.10(1) 


-2.5(1.0) 




diamond 


-3.0(2) 


0.10(2) 


-2.3(4) 




(7-expansion 


-3.92, -4.27 


0.126, 0.120 


-3.50, -4.12 




e-expansion 


-3.06, -4.99 


0.080, 0.134 


-3.14, -5.13 




impr-e-expansion 


-2.80, -3.64 


0.060, 0.089 


-2.86, -3.73 


2 


cubic 


-3.9(2) 


0.11(1) 


-3.5(1.0) 




diamond 


-4.1(4) 


0.10(2) 


-3.5(3) 




p-expansion 


-4.22, -4.54 


0.133, 0.128 


-3.85, -4.40 




e-expansion 


-3.39, -5.29 


0.089, 0.142 


-3.48, -5.44 


3 


cubic 


-4.1(1) 


0.11(2) 


-4.1(4) 




diamond 


-4.5(3) 


0.11(3) 


-4.0(4) 




^-expansion 


-4.29, -4.58 


0.134, 0.128 


-3.96, -4.45 




e-expansion 


-3.56, -4.55 


0.094, 0.144 


-3.66, -5.50 


4 


cubic 


-4.1(2) 


0.12(1) 


-4(1) 




diamond 


-4.7(2) 


0.10(2) 


-4.2(4) 




(jf-cxpansion 


-4.21, -4.46 


0.130, 0.125 


-3.92, -4.34 




e-cxpansion 


-3.64, -5.28 


0.096, 0.143 


-3.74, -5.43 




l/A'-cxpansion 


-11.12 


0.336 


-11.48 


8 


cubic 


-3.5(1) 


0.09(2) 


-3.8(5) 




diamond 


-4.0(1) 


0.05(5) 


-3.8(4) 




p-expansion 


-3.60, -3.72 


0.108, 0.103 


-3.44, -3.68 




e-expansion 


-3.48, -4.55 


0.093, 0.124 


-3.58, -4.68 




1 / Af-cxpansion 


-5.56 


0.118 


-5.74 


16 


cubic 


-2.4(1) 


0.06(1) 


-2.8(2) 




diamond 


-2.65(5) 


0.05(3) 


-2.7(3) 




5-expansion 


-2.46, -2.49 


0.072, 0.069 


-2.43, -2.52 




e-expansion 


-2.73, -3.19 


0.074, 0.088 


-2.81, -3.28 




1 / AT-expansion 


-2.78 


0.084 


-2.87 


32 


cubic 


-1.45(5) 


0.04(1) 


-1.8(3) 




diamond 


-1.50(5) 


0.04(1) 


-1.7(3) 




g-expansion 


-1.427, -1.429 


0.041, 0.040 


-1.45, -1.48 




e-cxpansion 


-1.73, -1.84 


0.047, 0.052 


-1.78, -1.90 




1 /A'-cxpansion 


-1.39 


0.042 


-1.43 


oc 
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TABLE VI. Three-dimensional 0{N) a model with nearest-neighbor interactions: lowest mo- 
ments of G{x) at iV = oo on the cubic and diamond lattice. 



moments 


cubic 


f.c.c. 


diamond 




i 


i 


3 



X 
m2 
Ml 

93,0 

m4 

94,0 
93,1 

me 

94,1 
96,0 

ms 

94,2 
96,1 



(1 + ^) 



12 
hz 



5040 I 1 I 1 \ 

I 10 840 j 



528 
51^ 



12 
77z 



2/3z 
6 







_3_ 

6z 



1+ ■ 



10 ^ 840 i 

-132 ("1 + ^^ 

5^ l,-^ ^ 44 ; 



Sz 
20 



-LIZ I Z 

2160 60480 



41184 



I ^ 286 ^ 3432 / 

(1 + tIo) 



\ 362880 /'i 



39 
" lB4z 



3z 
20 



37z"' 
7660 



10296 
6z^ 



240 
TIP 



I ^ 143 ^ 3432 / 







1110 
" 77z^ 



2/3z 



1 + ■ 



352 
5l^ 



27456 
"-553- 

23680 
231z^ 



llz , Szf. , z^ "\ 
60 946 10080 j 

(i + l§ + ^ + o(-^)) 



416 
"231z 



1 + 


3z 1 
20 


389z^ 
136080 


1 + 


868 ^ 


139z^ 
30888 




129z 
1480 


I 5z^ 
"T 444 



TABLE VIL Two-dimensional 0(N) a model with nearest-neighbor interactions: lowest mo- 
ments of G{x) at AT = 00 on the square, triangular, and honeycomb lattice. 



moments 



1712 

Ml 

*3,0 

m4 

94,0 

m6 

94,1 
t6,0 



square 



64 
7^ 



(1 + ^) 



(l+f + ^) 



40 
7^ 







triangular 



2 

3/3z 





(1 + ^) 



(l + f + ^) 



honeycomb 



Hi. 



3/3z 
4 



2 

64 
7^ 







" 4 64 4608)(-'- 
(l-#)(l+f)" 



62 



TABLE VIII. For various values of N, we report estimates of a obtained by our strong-coupling 
analysis, from the l/A'^-expansion, from the resummation of the g-expansion (see Section [1IID| ) (in 
this case we give two numbers corresponding to the two choices: resumming R{x) or i?(x)/x^), and 
from the O(e^) term of the e-expansion. In order to derive a from ai^, which is what is computed 
in the strong-coupling analysis, we used the following values of i^: ~ 0.59 for = 0; ~ 0.63 for 
AT = 1; ~ 0.67 for = 2; ~ 0.71 for A^ = 3; ~ 0.74 for A^ = 4; ~ 0.83 for A^ = 8; z^ ~ 0.91 
for A^ = 16; ~ 0.96 for A^ = 32. The errors we report for the strong-coupling estimates take into 
account all the analyses we performed. 



A'' s.c. -expansion 1/Af-cxpansion ly-expansion e-cxpansion 

"O 0.00(1) 0.0119, 0.0141 0.0109 

1 0.01(1) 0.0143, 0.0166 0.0130 

2 0.02(1) 0.0156, 0.0177 0.0140 

3 0.03(2) 0.0515 0.0160, 0.0179 0.0145 

4 0.03(2) 0.0386 0.0158, 0.0174 0.0147 
8 0.02(1) 0.0193 0.0139, 0.0148 0.0137 
16 0.009(3) 0.0096 0.0098, 0.0109 0.0109 
32 0.004(2) 0.0048 0.0058, 0.0059 0.0074 



TABLE IX. Estimates of erg, obtained by applying the CPRM to the series ^6,0 m2 on the 
cubic lattice. The errors reported in the Table take into account all the analyses we performed. 



N s.c. -expansion 1/Af-cxpansion e-expansion 

"O 0.01(1) 0.0134 

1 0.03(2) 0.0159 

2 0.04(2) 0.0171 

3 0.04(2) 0.0177 

4 0.036(10) 0.0491 0.0178 
8 0.024(4) 0.0245 0.0167 
16 0.013(2) 0.0123 0.0134 
32 0.0065(8) 0.0061 0.0091 
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TABLE X. Estimates of di from the analysis of the strong-couphng scries on the cubic and 
diamond lattice. The last column is our final estimate. An asterisk indicates that most of the 
approximants we considered are defective, or, in the cases where numbers are not shown, that 
all the approximants are defective, so that no estimate can be extracted. We also report results 



obtained by resumming the available terms of the p-expansion, from the 1/N calculation. 



N 




lO^'PA 


lO^DPA 


lO* lA 


iCdi 





cubic 
diamond 
g-expansion 


1.3(7) 
-0.9(5) 


* 


*0.4(9) 
*-1.2(2) 


1(1) 
-1.0(5) 
-1.31, -1.60 


1 


cubic 

(liaiiioiid 
(y-cxpaiibioii 


-1.6(8) 

-3(1) 


-3.1(7) 


*-1.7(l) 

*-.3.1(.'^) 


-1.7(5) 

-3(1) 

-l.T)!), -1.89 


2 


cubic 
diamond 
g-expansion 


-2.2(3) 
-3(1) 


=1= 

-2.8(3) 


-2.3(1) 
* -3.7(9) 


-2.3(2) 
-3(1) 

-1.72, -2.01 


3 


cubic 

diamond 
(/-expansion 


-2.4(3) 
-2.2(9) 


* 

-2.6(3) 


-2.5(1) 
-2.6(7) 


-2.5(2) 

-2.6(3) 
-1.77, -2.03 


4 


cubic 
diamond 
g-cxpansion 
1 / Af-expansion 


-2.4(3) 
-5(3) 


* 

-2.3(3) 


-2.5(2) 
*-4(3) 


-2.5(2) 
-2.5(5) 
-1.76, -1.99 
-5.12 


8 


cubic 
diamond 


•-2.0(4) 
-5(3) 


* 

-2(2) 


-2.1(2) 
*-4(l) 


-2.1(2) 
-3(2) 




g-expansion 
1 /AT-expansion 








-1.55, -1.68 
-2.56 


16 


cubic 
diamond 
g-expansion 
1 /AT-expansion 


-1.3(3) 
-3(3) 


* 

-1.2(8) 


-1.4(2) 
* 


-1.4(2) 
-1.2(8) 
-1.10, -1.15 
-1.28 


32 


cubic 
diamond 
g-expansion 
1 /AT-expansion 


-0.7(3) 
-1(1) 


* 

-0.4(2) 


-0.7(3) 
* 


-0.7(2) 
-0.5(3) 
-0.66, -0.67 

-0.64 


oo 














64 



itr 

I 



1 

2 
3 
4 
5 
6 
6 
7 
8 
8 
8 
9 
9 
9 
10 
10 
10 

11 

11 
11 
11 

12 
12 
12 
12 
12 
12 
13 
13 
13 
13 
14 
14 
14 
14 
14 
14 
14 
15 
15 
15 
15 
15 
15 
15 
15 



KI. Coefficients of the strong-coupling expansion of G ^{x) on the cubic lattice. The 
re of each equivalence class is chosen by xi > X2 > > 0. / indicates the order. 



N = 



N =1 



N : 



N : 



N : 



AT = 16 



-1 



-1 



-1 



30 



26 



24 



114 

5 



22 



68 

3 



-13 



122 
15 



6 



158 
" 35 



11 

' 3 



49 
" 45 



366 
2 



4204 
15 





479 

2 
_ 1 
2 



175 
_ 16 
25 



602 

3 
_ 2 
3_ 



135 
— 10 
27 



-197 



33604 
316 



1123 
16 



44812 
875 



1793 
45 



132961 
13366 



5022 

24 
1 



106 

-16 



4 
65 



-1 



875 
18576 



25 



4455 
_ 1912 
243 



-2889 
6 

-1 







240 
-57 
4 
_ 1 

4 



9626 
408 

25 
56 
125 



135 
50 

" 3 

_5 
9 



2940271 
40096 
778 
81 
115 
243 



76062 

258 
116 



211434604 
4725 
962 
3 

-24 



3833513 
120 
24851 
72 
683 
12 



8484905796 
336875 
1878984 
6125 
45272 



21078 

35786 
135 
1340 
27 



1969903274 
200475 
10005158 
120285 
49060 



-45357 
72 

-15 




2874597004 
166926 
-240 

-24 



4320 
2635 

8 
133 

3 

_ 1 
2 



949030894596 
153278125 
7370064 

21875 
1142304 
21875 
96 

125 



2835 
8636 

27 
4369 

81 



34414485049 
70366725 
915928 

6561 
605717 
19683 
160 
243 



1230462 
2460 
1944 









101355262012 
165925 
264288 
45 

-896 
-24 



77441167 
180 
470691 

90 
43417 
36 
393 



49050932335452 
153278125 
49887963408 
11790626 
36397888 
30626 
39096 
625 
1104 
625 
192 

625 



26831414 
106 
1379716 

405 
423976 
405 
616 
9 

_64 
27 

_4 
9 



2184687177202 
23456575 
3960090404 

5412826 
427970066 
1082666 
35456 
729 
6600 
2187 
400 

729 



-745189 
-678 
-476 




1629512844964 
6081075 
28728 
5 

-1120 

-24 



2652479497 
20160 
614743 

96 
437317 
288 
1337 
24 



11642408503972 
153278125 
903833736 

153125 
241988968 
153125 
15(10984 
21875 



945 
704326 

135 
616348 

405 

5864 
81 



470143113389 
164189025 
600662842 

360855 
765424036 
1082565 
868672 
19683 



20787102 
17378 
29088 


32 
-2 




140597688722408 
14189176 
18013984 
189 

-20744 



-1456 
-144 
-32 



245838046393 
40320 
1330017637 
17280 
8258723 
360 
54377 

24 
2425 

8 
206 
3 

_ 1 
4 



298587843693288 
69671875 
43519901193056 

766390625 
1192015283352 
58953125 
278755152 

109375 
42568048 
109375 
1427216 
15625 
288 



46330187578 
14175 
1801726602 
42626 
21049024 
1215 
628336 

243 
103640 
243 
8294 
81 
_16 
27 



15116615635191578 
16264713476 
36638833424418 
6966305775 
253843671008 



48715425 
253847360 

177147 
52536352 

177147 
4477778 



59049 
1280 
21-^7 



-12672757 
-48624 
-16428 
30 
-400 
-10 
-1 




2597638257008408 
638512875 
7317712 

63 
164436 
6 

-232 
-1752 
-48 







3549389785799 
1935360 
317601607 

2880 
605816971 
17280 
8199 
16 
192203 
72 
1681 
16 
_ 13 
8 

_ J_ 
16 



64161784918165784 
65143203125 
655516678768 

6016626 
1403163936156 
42109376 
2141976 
3125 
11182524776 
3828125 
90832 

625 
11936 
3125 
648 
3126 



25033873577 
42526 
92268604 

1216 
110231638 
3646 
20964 
27 

10592632 
3645 
506 

3 
151 
27 
— 28 
81 



3051121951)55173211 
24138249510375 
2363546376656 

146146276 
4917927730652 
438438826 
4066162 
6661 
133926221684 
87687765 
953338 
6561 
50545 



6561 
11840 
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